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Jean  Rostand 


PREFACE 


The  study  of  detonation  waves  has  been  the  object,  for 
twenty  years  or  so,  of  numerous  theoretical  as  well  as  experi¬ 
mental  investigations. 

Of  all  the  as  yet  unsettled  problems  in  this  vast  field, 
we  have  examined  the  initiation  of  detonation  by  shock  wave  in 
a  solid  granular  explosive.  Whereas  at  the  time  when  this 
study  was  begun  the  behavior  of  homogeneous  explosives  was 
quite  well  explained,  there  had  been  only  a  few  experimental 
studies  dealing  with  the  behavior  of  granular  explosives.  In 
order  to  express  these  results  — -  which  are  sharply  different 
from  those  furnished  by  the  homogeneous  explosives  —  theoret¬ 
ically,  we  have  worked  out  a  mathematical  scheme  of  the  mechan¬ 
ism  of  Initiation  of  the  detonation  by  shock  wave. 

I  would  like,  on  the  occasion  of  editing  this  work,  to 
pay  homage  to  the  memory  of  Professor  Ribaud,  who  died  last 
autumn,  and  to  thank  Professor  Rocard  for  agreeing  to  devote 
his  attention  to  this  quite  special  problem  of  fluid  mechanics, 
despite  his  numerous  responsibilities.  I  also  thank  M.  Robert, 
Director  of  the  Atomic  Energy  Commission,  for  authorizing  me 
to  carry  out  this  work  in  his  laboratories,  as  well  as  y.  Bar- 
gulllet,  Director  of  the  Research  Center  at  Vajours,  for  all 
the  facilities  which  he  has  placed  at  my  disposal. 

I  want  to  express  my  profound  gratitude  to  M.  Berger, 
Chief  Explosives  Engineer,  for  his  enlightened  comments,  his 
constructive  criticism  and  the  equipment  which  he  was  so  kind 
to  make  available  to  me. 

Finally,  I  must  express  my  thanks  to  M.  VId.art  and  Mme 
Prouteau  for  the  programming  and  treatment  of  the  mathematical 
scheme  on  the  computer,  and  to  M.  Lezaud  and  his  team  for  their 
constant  devotion  during  the  preparation  and  execution  of  the 
experiments . 
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INTRODUCTION 


Historically,  it  seems  that  the  initiation  of  detonation 
by  Impact  in  an  explosive  wa3  discovered  accidentally.  It  was 
probably  the  accidents  which  were  unexplainable  at  that  time 
that  drew  the  attention  of  researchers  to  this  phenomenon. 

At  the  present  time  the  procedure  is  quite  frequently 
employed  by  laboratories  specializing  in  the  study  of  explosives, 
and  also  of  certain  detonators  where  the  effect  of  impact  is 
added  to  the  thermal  and  shock  phenomena.  To  this  end,  the 
impact  is  brought  about  by  a  metallic  plate  generally  projected 
by  means  of  a  sheet  of  explosive. 

Though  frequently  used,  the  phenomenon  itself  has  been 
studied  only  to  a  limited  extent.  One  hae  been  content  with 
verifying  that  it  was  indeed  the  Impact  which  brings  about 
the  detonation,  and  with  roughly  determining  the  characteris¬ 
tics  to  be  used:  the  velocity  thresnold  for  a  given  projectile 
and  explosive.  Nevertheless  it  should  be  pointed  out  that  these 
studies  had  an  essentially  practical  objective  —  testing  of 
safety  —  and  were  not  oriented  toward  basic  research. 

It  was  only  in  1961  (Whitbread,  Drown*)  that  certain 
parameters  specific  to  this  type  of  initiation  of  detonation 
were  stated  precisely. 


*S.  M.  Brown,  E.  b.  Whl thread .  Les  Ondes  de  p£tnnatlon 
(Detonation  Waves),  1962,  Publications  of  the  CNR3  (Centre 
National  de  Recherche  Scientlf lque ;  National  Center  for  Scien¬ 
tific  Research),  No  109. 
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While  the  problem  of  detonation  produced  by  impact  has 
received  only  limited  experimental  treatment,  we  have  numerous 
results  relating  to  a  very  similar  phenomenon,  thau  of  detona¬ 
tion  produced  by  shock  (with  barrier).  In  Doth  cases  it  is 
the  Induced  shock  wave  which  leads  to  the  detonation  of  the 
explosive. 

These  two  studies  have  many  points  in  common  with  the 
phenomenon  of  transition  from  deflagration  to  detonation. 

During  detonation  produced  by  impact,  the  physical 
phenomena  involved  seem  to  differ  according  to  the  intensity 
of  this  impact.  In  the  case  of  a  very  weak  induced  shock, 
e.g.  that  generated  by  Impact  in  sensitivity  tests  of  explosives 
(fall  of  a  mass  weighing  a  few  kilos  from  a  height  of  approx¬ 
imately  one  meter),  the  detonation  seems  to  be  of  thermal 
origin.  In  particular,  the  “hot  spots"  hypothesis  seems  t  > 
have  received  a  clear  experimental  confirmation  by  the  wo  s.  of 
Bowden  and  his  collaborators.  Likewise,  when  the  detonation 
is  brought  about  by  heat  or  by  light,  the  explosive  is  subjected 
to  a  physical-chemical  evolution  in  which  the  thermal  phenomena 
seem  to  play  a  large  role. 

This  is  not  the  case  when  a  shock  wave  is  generated. 

Taken  as  a  whole,  the  phenomenon  is  so  fast  (a  few  microseconds) 
that  in  order  to  explain  the  behavior  of  the  explosive,  one 
must  call  upon  considerations  other  than  those  of  thermal  con¬ 
duction  and  convection. 

In  the  present  study  we  shall  deal  only  with  the  behavior 
of  explosives  subjected  to  high- intensity  shocks  (several  kilo- 
bars  ) . 


To  bring  about  such  shocks,  various  devices  may  be 
employed.  The  most  common  shock  generators  are  the  explosives 
themselves.  However,  in  that  case,  the  shock  is  a  reactive  one 
and  in  order  to  transform  it  into  inert  shock,  it  Is  necessary 
to  proceed  through  the  intermediary  of  a  barrier.  This  pro¬ 
cedure  has  been  tried  out  in  numerous  laboratories.  It  has 
made  it  possible  to  obtain  valuable  data  regarding  this 
phenomenon  of  detonation  initiation,  but  has  the  serious  defect 
that  it  Introduces  a  pressure  signal  Into  the  explosive  which 
is  not  too  well  known.  Hence  the  interpretation  of  the  results 
is  a  difficult  undertaking. 

The  other  method  of  initiation  by  shock  wave  is  the 
Impact  of  a  projectile  launched  at  a  high  speed,  the  explosive 
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serving  as  target.  This  procedure  presents  the  advantage  of 
assuring  a  shock  of  constant  characteristics  during  a  certain 
period  of  time.  However,  few  experimental  results  have  been 
published . 

The  explosives  subjected  to  shock  do  not  behave  In  an 
Identical  manner.  The  homogeneous  explosives  (liquids,  single 
crystals)  exhibit  a  period  of  "incubation"  prior  to  the  Initi¬ 
ation  of  the  detonation,  while  the  heterogeneous  explosives 
(granular  solids)  do  not  exhibit  such  a  period.  The  detailed 
examination  of  these  differences  Is  completed  by  a  study  of 
the  various  physical  mechanisms  which  must  be  called  upon  to 
explain  the  initiation  of  the  detonation;  this  forms  the  object 
of  Chapter  2. 

The  next  chapter  is  devoted  to  the  experimental  apparatus 
used  by  us,  both  for  the  realization  of  a  correct  Impact  of 
known  characteristics,  and  for  the  study  of  the  shock- 
detonation  transition  In  the  xploslve  serving  as  target. 

These  devices  have  the  peculiarity  of  assuring  an  almost 
one-dimenslonal  phenomenon  in  the  entire  zone  of  measurement. 
This  makes  It  possible  to  obtain  results  that  may  oe  directly 
compared  with  those  furnished  by  calculation. 

Next,  the  theoretical  aspect  of  the  problem  Is  approached 
Strictly  speaking,  If  we  want  to  take  into  account  that  the 
explosive  Is  not  a  homogeneous  body  but  Is  made  up  of  essentiall 
epnerical  granules  In  contact  with  each  other,  then  the  hydro- 
dynamic  phenomenon  alone  Is  nonstationary  tridimensional. 
Noreover,  the  medium  Is  reactive.  Hence  It  is  necessary  to 
make  use  of  chemical  kinetics  and  the  various  equations  of 
state  of  the  constituent  substances  (solids)  or  of  those  formed 
In  the  course  of  the  reaction  (generally  gases). 

Hence  a  simplifying  scheme  Is  necessary.  Its  only  pur¬ 
pose  Is  to  account  fully  for  the  overall  phenomenon. 

.ihile  the  presently  available  schemes  do  give  a  good 
explanation  of  the  behavior  of  homogeneous  explosives,  they  do 
not  make  It  possible  to  explain  the  characteristic  behavior  of 
granular  explosives.  In  our  opinion  this  seems  to  be  due  to 
two  reasons:  first,  the  solid  and  gaseous  phases  cannot  be  at 
the  same  temperature  for  the  same  pressure;  and  second,  the 
Arrhenius- type  law  of  chemical  kinetics  is  much  too  violent  to 
initiate  the  reaction.  A  mathematical  model  which  takes  Into 
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account  these  two  remarks  —  differentiation  of  the  temperatures 
of  the  solid  and  gaseous  phases  during  the  reaction;  a  law  of 
chemical  kinetics  that  is  more  progressive  than  Arrhenius' 
exponential  law  —  was  constructed.  It  furnishes  a  rather 
satisfactory  agreement  with  the  various  known  experimental 
results. 
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Chapter  1 

STUDY  OF  THE  SHOCK  WAVE  RESPONSIBLE  FOR  THE  INITIATION 

OF  DETONATION 


1.1.  RECAPITULATION  OF  THE  HYDRODYNAMIC  EQUATIONS  OF  SHOCK. 

APPLICATION  TO  SOLIDS. 

Let  us  consider  a  fluid  initially  in  state  (0)  traversed 
by  a  shock  wave  of  absolute  velocity  U,  downstream  of  which  the 
state  of  the  fluid  is  denoted  by  (1). 

In  a  one-dlrrenaional  flow,  the  equation  of  the  conser¬ 
vation  of  mass  is  written  as 

pjv  -  «.)  •  P»(U  -  u,) 

The  momentum  theory  furnishes  the  equation 


ft  -  P.  •  -  u.)  (u,  -  u.) 

and  the  conservation  of  energy: 


Pi  ui  *  p. 


p.<U  -  «.)  [e,  -  E.  ♦  |(u»  .  u»)J 


In  the  following  applications,  p0  is  negligible  compared 
with  p,  and  u0  may  be  taken  as  zero. 


After  rearrangement,  the  above  system  of  equations  may 
be  written  in  the  following  form  which  can  be  used  more 
directly: 

p.  u  •  Pl(u  -  u  ) 

p  •  p,  u  u 


El  -  E. 


6 


For  a  variable  shock  intensity,  we  obtain  in  the  p,  v 

plane  (v  =  -)  the  dynamic  adiabatic,  and  in  the  p,  u  plane  the 

P 

shock  polar  of  the  body,  while  the  equations  of  state  and  of 
Internal  energy  furnish  the  two  complementary  equations  neces¬ 
sary  for  the  complete  determination  of  the  various  variables  of 
the  shock  as  a  function  of  one  of  them  chosen  aB  parameter. 

The  equation  of  state  of  a  solid  may  be  written  as: 

p  •  p4  ♦  c  tt*  -  T.) 

In  this  way  the  equation  of  isothermal  compressibility 
is  stated  in  a  clear  manner. 

More  precisely,  we  have  chosen  for  all  subsequent  cal¬ 
culations  the  term  given  by  the  equation  of  Pack,  Evanj  and 
James  [14]: 

p‘  ‘  *  ("t J*  fxp  [‘  *  (7/]  ]  -  *] 

and  taken  the  coefficient  g  =  as  constant. 

OT  v 

With  this  equation  of  state,  an  isentroplc  defined  by 

dE  ♦  p  dv  •  0 


with 

dE  -  c.  dT  ♦  (T  -  p)  dr 

where  cY  is  assumed  to  be  constant,  leads  to 

P  *  P,  -  f  T.  ♦  f  T,  (tp  J  J-  <v,  .  y)J 

with  subscript  1  indicating  the  pole  of  this  isentroplc. 
The  temperature  is  given  by 


T  •  T,  up  J  J-  -  v)  J 

and  the  velocity  of  sound,  defined  as  follows: 
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is  then  given  by  the  expression: 


'  •  t*(^/  |  [»p  *(•■(£)/][•♦>  %t]  «4rT*  “"O'  - «’] 

The  equations  of  the  shock  furnish  the  dynamic  adiabatic 

p(v): 


p 


-  (  T.(v  -  v.)  ♦  /* 

_  M 


v) 


p4  dv 


This  has,  for  pole,  the  point  p  =  0,  v  =  vc  of  Clapeyron' 
diagram. 

The  velocity  of  propagation  of  the  shock,  U,  and  the 
material  velocity,  u,  are  then  given  by 


u* 


The  temperature  is  Immediately  deduced  from  p,  by  means  of  the 
equation 


X  •  T.  *  P  -I-£i 

In  the  subsequent  numerical  applications,  we  have  used, 
for  the  sake  of  simplifying  the  calculations,  linear  relation¬ 
ships  between  the  various  parameters. 

In  particular,  we  know  ([3],  p  261)  that  there  exists, 
in  a  very  good  approximation  and  in  a  relatively  large  zone,  a 
linear  relationship  between  the  velocity  of  the  shock  (U)  and 
the  material  velocity  (u):  U  =  A  +  Bu. 

A  second  approximation  was  employed  for  determining  the 
velocity  of  sound  as  a  function  of  the  ratio  p/p0, 


This  last  law,  which  is  the  linearization  of  the  pre¬ 
ceding  theoretical  results,  has  been  experimentally  proved  in 
a  quite  satisfactory  manner  i2]. 

These  two  simplifications  have  permitted  us  to  reduce 
considerably  the  time  of  calculation  and  to  preserve  the  pos¬ 
sibility  of  varying  the  different  parameters  while  insuring  a 
satisfactory  accuracy. 

1.2.  PRESENTATION  OF  THE  DIFFERENT  MODES  OF  INITIATION  OF  THE 
SHOCK 

The  procedure  most  commonly  used  for  obtaining  shocks 
of  high  intensity  employs  an  explosive.  This  explosive  Induces 
in  a  barrier  an  inert  shock  which  is  subsequently  transmitted 
to  the  receiving  explosive.  Two  procedures  of  initiation  may 
be  considered:  frontal  shock  and  lateral  shock. 


e  Ik 


I.  © 


$ - 

r#c.p«.ur 


lap act 


(1)  Generating  explosive;  (2)  Receptor;  (3)  Frontal 
shock;  (A)  Barrier;  (5)  Lateral  shock;  (6)  Gener¬ 
ator;  (7)  Target. 

For  obvious  reasons  of  convenience  of  interpretation, 
only  the  frontal  shock  is  employed  in  these  studies.  So.i.e  tests 
have,  indeed,  been  carried  out  by  lateral  shock  [6],  but  to  our 
knowledge  no  theoretical  study  has  been  subsequently  undertaken 
of  the  trfo-dimensional  nonstationary  character  of  the  phenomenon 
(three  variables:  two  of  space  and  one  of  time).  By  contrast, 
the  frontal  shock  may  be  one-dimensional  nonstationary  (only 
two  variables:  one  of  space,  the  other  of  time)  If  one  adheres 
to  certain  experimental  conditions.  Hence  Its  physical  inter¬ 
pretation  and  Its  being  put  into  an  equation  form  are  a  priori 
easier. 

Another  procedure  for  the  Initiation  of  shock  consists 
in  the  Impact  of  a  projectile  on  a  target.  This  method  has 
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made  it  possible  to  bring  about  the  most  intense  shocks  up  to 
now  (up  to  10  megabars)  Li]  il  1 3 3 •  If  certain  experimental  pre¬ 
cautions  similar  to  those  employed  in  the  case  of  frontal  shock 
are  taken,  it  can  lead  to  a  one-dimensional  nonstationary 
phenomenon. 

1.3.  DETERMINATION  OF  THE  PARAMETERS  OF  THE  SHOCK  INDUCED  IN 
AN  INERT  RECEPTOR 

After  recalling  the  methods  used  in  this  determination 
when  the  tested  body  is  an  inert  solid,  we  shall  discuss  the 
validity  of  these  procedures  when  the  receptor  is  an  explosive 
substance. 

1.3.1.  Experimental  Methods 

The  intensity  of  the  shock  induced  in  the  receptor  is 
determined  by  the  equality  of  the  pressures  (p)  and  of  the 
material  velocities  (u)  on  both  sides  of  the  interface  (barrier- 
receptor  in  the  case  of  frontal  shock;  projectile-target  in  the 
case  of  Impact).  Hence  it  is  advantageous  to  employ  the  (p,  u) 
plane  for  graphical  representation.  In  this  plane  the  shock 
polar  p(u)  of  the  barrier  (or  the  projectile)  is  assumed  known. 
Furthermore,  it  may  be  determined  previously  by  means  of  analo¬ 
gous  experiments,  provided  that  the  receptor  is  made  of  the 
same  material  as  the  barrier  (or  the  projectile) 

Experimentally  it  Is  sufficient  to  determine  two  param¬ 
eters  of  the  shock.  The  hydrodynamic  equations  permit  calcu¬ 
lating  the  other  variables.  The  speed  of  propagation  of  the 
shock  is  determined  by  measuring  the  time  elapsed  during  the 
passage  of  the  shock  between  two  points  situated  at  a  known 
distance  from  each  other.  The  other  parameter  measured  is  the 
free  surface  velocity  ug^  (velocity  acquired  by  the  maes  after 

expansion  to  atmospheric  pressure).  The  latter,  in  a  good 
approximation,  is  equal  to  twice  the  material  velocity  L 3 » 

P  263 J . 


Various  methods  are  employed  for  carrying  out  these 
measurements : 

The  optical  method  of  chambers  filled  with  argon  (or 
xenon)  l  3 »  P  267]  employs  the  Intense  luminosity  of  this  gas 
under  impact.  By  means  of  the  slit  camera  the  instant  of 
closure  of  these  chambers  is  recorded. 

The  method  of  electric  probes  is  b^sed  on  the  short  cir¬ 
cuit  which  may  be  obtained  by  the  movement  of  the  material  as 
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soon  as  the  shock  passes.  The  moment  of  closing  of  this  cir¬ 
cuit  breaker  placed  into  an  electronic  circuit  is  registered  by 
means  of  a.  oscilloscope. 

The  wire-  (or  luminous  slits)  method  consists  in  measur¬ 
ing  the  rate  of  displacement  of  a  reflecting  surface  by  means 
of  the  apparent  rate  of  displacement  of  the  image  of  a  fixed 
wire  (or  a  luminous  slit)  [5]  [153* 


*1 

•1 


Intarfac*  ?0  -0 

“o  "0 


conditions  Initial** 
(avant  I'arrivi*  k  ® 
I' intarfac*) 


Interface 


(1)  Shock;  (2)  Shock  polar  of  the  barrier;  (3)  Expansion; 

(4)  Initial  conditions  (before  arrival  at  the  interface); 

(5)  Case  of  shock;  (6)  Case  of  expansion. 


1.3.2.  Case  of  the  Frontal  Shock 

The  measurement  of  the  free  surface  velocity  of  the 
barrier  makes  It  possible  to  place,  In  plane  (p,  u),  the  repre¬ 
sentative  point  (p1t  Uj)  on  the  shock  polar  of  this  body,  while 

the  measurement  of  the  shock  velocity  U  in  the  receptor  medium 
permits  drawing  the  straight  line  having  the  slope  p0  U  (where 

(>0  1  s  the  initial  specific  mass  of  the  receptor).  Depending 

on  the  position  of  this  straight  line  with  respect  to  point 
(Pj*  Uj),  a  shock  or  a  beam  of  expansion  waves  rises  In  the 

barrier,  while  a  shock  is  propagated  in  the  receptor. 

This  method,  known  as  the  "graphic"  method,  has  been 
used  for  determining  the  dynamic  adiabatic  of  the  majority  of 
common  metals. 

It  should  be  noted  that,  strictly  speaking,  the  point 
representative  of  the  shock  (P2»  rises  In  the  barrier 

Is  situated  on  the  curve  symmetrical  with  respect  to  straight 
line  u  =  u,  of  the  shock  polar  having  Pj ,  u^  for  pole,  and  not 

on  the  symmetric  of  the  shock  polar  having  p  =  0,  u  =  0  for 
pole.  The  case  is  the  same  when  we  go  from  (pj»  u1 )  to  (pg, 

U2)d  via  an  expansion;  strictly  speaking,  the  latter  is  lsen- 

troplc.  However,  in  the  case  of  solids  it  may  be  estimated 
that  these  curves  coincide  in  a  large  pressure  range. 

1.3.3.  Case  of  Impact 

The  measurements  carried  out  are:  the  impact  velocity 
of  the  projectile  and  the  shock  velocity  in  the  receptor. 

The  conditions  of  equal  pressure  and  mass  velocity  at 
the  interface  are  given,  in  plane  (p,  u),  by  the  intersection 
of  the  shock  polar  of  the  projectile  and  the  straight  line 
having  the  slope  pn  U  (U  =  the  shock  velocity  in  the  receptor 

of  initial  specific  ma3s,o0). 

In  contrast  to  the  frontal  shock  method,  the  point 
representative  of  the  shock  in  the  receptor  is  situated 
strictly  on  the  curve  which  is  symmetrical  with  respect  to  the 
straight  line  u  =  u1/2  of  the  shock  polar  of  the  projectile  of 
pole  p  r-  0 ,  u  =  0 ) . 


12 


(1)  Symmetric  of  the  shock  polar  of  the 
projectile;  (2)  Impact  velocity. 

1.3.4.  Shape  of  Pressure  Signal  Induced  in  the  Receptor 

The  shock  Induced  in  the  barrier,  and  consequently  in 
the  receptor,  is  not  strictly  constant.  In  effect,  in  addition 
to  the  uncertainty  with  respect  to  the  reaction  zone  and  the 
shape  of  the  pressure  signal  in  this  zone,  the  expansion  of 
detonation  products  constantly  perturbs  the  Induced  shock. 

Hence  the  intensity,  shape  and  decay  in  time  of  the  pressure 
signal  Induced  in  the  receptor  cannot  at  the  present  time  be 
determined  accurately:  they  depend,  for  the  same  generating 
explosive,  on  the  nature  and  thickness  of  the  barrier. 

Hence  this  generator  is  not  too  satisfactory  since  by 
its  very  design  it  makes  one  feel  certain  that  it  Induces  in 
the  receptor  an  essentially  variable  pressure  signal;  this 
signal  decays  in  space  and  time  Independently  of  any  dissipat¬ 
ing  phenomenon.  Hence  it  is  indispensable  to  calibrate  the 
generator  by  studying  the  decay  of  the  shock  in  space  and  time 
as  a  function  of  the  barrier  thickness;  then  it  is  possible  to 
determine  the  initial  intensity  of  the  shock  for  a  given 
barrier  thickness,  roughly  estimate  the  shape  of  the  pressure 
signal  and  represent  it,  at  a  given  Instant,  by  an  exponential 
law  or  by  a  linear  law  (signal  of  triangular  shape). 

In  the  case  of  impact  the  phenomenon  is  different.  In 
effect,  the  shocks  induced  in  the  projectile  and  in  the  target 
are  of  constant  intensity  except  for  the  effect  of  any  dissi¬ 
pation  phenomena,  as  long  as  a  free  surface  has  not  been 
attained  by  one  of  them.  Hence,  since  the  receptor  may  be 
chosen  to  be  of  a  sufficient  length  to  rule  out  this  incon¬ 
venience,  it  is  the  thickness  of  the  projectile  which  determines 
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the  time  during  which  the  shock  Induced  In  the  explosive 
remains  constant. 


This  nevertheless  supposes  that  at  the  moment  of  impact 
the  pressure  is  zero  and  the  velocity  Is  equal  to  a  constant 
at  all  points  of  the  projectile. 


For**  du  signal  da  praaalonQ) 
dana  la  barritra,  au  covra 
du  taapa. 


Cat  du  choc  induit  avac 


barntra  & 


Foraa  du  signal  da  prassion  dana  (3) 
la  cibla,  au  coara  du  taapa 


Caa  d«  1' impact  (JO 

(1)  Generating  explosive;  (2)  Barrier;  (3)  Receptor;  (4)  Prog¬ 
ress  diagram;  (5)  Case  of  shock  induced  with  barrier;  (6)  Decay 

(7)  Form  of  pressure  signal  in  barrier  In  the  course  of  time; 

(8)  Target  explosive;  (9)  Form  of  pressure  signal  In  target  In 
the  course  nr  time;  (10)  Case  of  Impact. 


Hence  this  mode  of  shock  generation  presents  two  im¬ 
portant  advantages  compared  with  frontal  shock:  the  shock 
Induced  in  the  receptor  is  initially  uniform,  and  the  time 
during  which  this  signal  is  permanent  depends  on  the  thickness 
of  the  projectile. 

1.3.5.  Generators  of  Shock  by  Impact 

Now  we  are  in  a  position  to  discuss  the  choice  of  the 
various  experimental  devices  for  realizing  the  impact. 

The  simplest  solution  consists  in  employing  a  cannon 
whose  projectile  exhibits  a  plane  face  in  the  front  so  as  to 
achieve  a  correct  Impact  on  a  sufficiently  large  surface  area. 
The  phenomenon  brought  about  in  the  target  is  thus  one¬ 
dimensional  if  care  is  taken  that  the  expansions  originating 
from  the  edges  of  the  generator  or  the  receptor  do  not  inter¬ 
vene  in  the  zone  where  the  measurements  are  carried  out. 

Nevertheless  the  classic  projectiles  have  a  rather 
complicated  shape  and  it  is  difficult  to  define  their  "useful" 
thickness.  On  the  other  hand,  their  use  does  not  permit  the 
variation  of  this  thickness  in  a  simple  manner.  The  use  of 
guide  shoes  considerably  complicates  the  mounting  and  makes  it 
necessary  that  one  make  sure  that  they  do  not  modify  the 
phenomenon  after  Impact.  This  method  has  been  successfully 
used  by  Brown  and  Whitbread,  enabling  them  to  obtain  the  only 
results  specific  to  this  mode  of  generation  of  detonation 
published  to  date,  even  though  the  one-dimensional  character 
of  the  phenomenon  was  probably  not  respected  [4]. 

It  nevertheless  seems  preferable  to  look  for  a  method 
of  projection  which  would  make  it  possible  to  choose  from 
among  a  large  variety  of  projectiles,  impact  velocities  (these 
two  parameters  determining  the  intensity  of  the  induced  shock) 
and  projectile  thicknesses  (which  defines  the  duration  of 
application  of  the  shock). 

The  procedure  known  as  the  "plate  lifting"  meets  these 
requirements  rather  well.  It  consists  of  laterally  projecting 
a  thin  metallic  plate  (a  few  millimeters  thick)  by  means  of  an 
explosive  whose  detonation  is  initiated  in  such  a  way  that  it 
propagates  parallel  to  the  plane  of  the  plate.  It  is  noted 
experimentally  that  the  latter,  after  a  zone  of  acceleration, 
becomes  essentially  plane  and  assumes  a  position  which  makes  a 
constant  angle  <p  with  its  original  direction.  Moreover,  this 
device  makes  it  possible  to  select  for  the  projectile  the 
three  fundamental  parameters  which  determine  the  shock  Induced 
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In  the  target.  In  view  of  the  diversity  of  possibilities  which 
it  offers,  we  have  preferred  this  device  to  all  others. 

Another  method  of  projection  of  the  plate  is  the  frontal 
method,  and  the  velocities  realized  in  this  manner  are  greater 
than  those  furnished  by  the  lateral  method.  In  addition  to 
presenting  the  same  advantages  as  the  preceding  method,  it 
makes  it  possible  to  explore  a  higher  velocity  range.  he  shall 
see  that  this  advantage  is  hard  to  make  use  of  in  this  study, 
since  then  the  steady  state  of  detonation  is  attained  too  fast 
in  the  target. 


Nttfcod*  l«*4r«1*  (|)  Wthod*  fr0Bt,i, 

(r»Uv«Mnt  da  alaqaa)  (4) 


(1)  Explosive;  (2)  Lifted  metal  plate;  (4)  Lateral 

method  (lifting  of  plate);  (4)  Frontal  method; 

(5)  Metal  plate. 

It  should  be  noted,  however,  that  the  stress  undergone 
by  the  metal  during  Its  being  set  In  motion,  especially  In  the 
device  for  the  lifting  of  the  plate,  may  Influence  its  subse¬ 
quent  behavior  by  allowing  the  survival,  within  the  metal,  of 
pressure  and  velocity  gradients.  However  after  a  trajectory  of 
a  few  centimeters  In  the  case  of  a  plate  having  a  thickness  of 
the  order  of  one  mm  or  a  few  mm,  it  can  be  assumed  that  these 
gradients  are  too  small  to  be  able  to  modify  the  behavior  of 
the  plate  upon  impact.  The  theoretical  study  of  the  following 
chapter  will  make  it  possible  to  affirm  this  hypothesis  in  a 
more  concrete  manner. 

1.4.  APPLICATION'  TO  THE  CASE  V.HERE  THE  RECEPTOR  IS  AN  EXPLOSIVE 

The  preceding  study  on  the  determination  of  the  param¬ 
eters  of  the  shock  Induced  in  a  receptor  implies  two  hypotheses 

--  The  media  in  contact  with  each  other  are  homogeneous; 
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—  No  process  other  than  that  of  hydrodynamics  perturbs 
the  phenomenon. 

Let  us  now  see  how  valid  these  hypotheses  are  In  the 
case  where  the  receptor  Is  an  explosive. 

I.4.1.  Are  the  Media  In  Contact  Homogeneous? 

This  hypothesis  Is  not  verified  when  the  receptor  Is  a 
solid  body  constituted  of  grains  of  variable  size  and  of  differ¬ 
ent  chemical  composition  having  a  greater  or  lesser  number  of 
gaseous  Inclusions,  depending  on  the  density  of  the  charge. 

The  behavior  of  such  a  mixture  under  shock  cannot  be  calculated 
theoretically.  Hence  It  Is  questionable  whether  the  hydro- 
dynamic  equations  of  shock  are  applicable. 

Then  the  problem  arises  as  to  whether  It  Is  legitimate 
to  speak,  In  the  case  of  this  mixture  of  substances,  of  a 
dynamic  adiabatic  when  the  experimentally  determined  values 
(shock  velocity,  free  surface  velocity,  and  even  charge  density) 
are  but  averages  of  the  values  of  each  of  the  constituents. 

We  know,  for  example,  that  when  the  charge  density  is  low,  the 
dynamic  adiabatic  which  Is  determined  experimentally  Is  prac¬ 
tically  that  of  the  gas  contained  In  the  explosive. 

Recent  studies  on  porous  substances  l10][11]  have 
revealed  notable  deviations  of  behavior  from  that  of  homogeneous 
solids . 


This  heterogeneous  character  of  the  receptor  confers  on 
the  front  of  the  wave  which  propagates  In  the  receptor  a 
sinusoid  form  In  perpetual  evolution.  Its  propagation  is  no 
longer  one-dimensional  and  related  phenomena  (for  example  wave 
convergence)  may  become  dominant  during  the  initiation  of  the 
detonation . 

1.4.2.  Is  the  Phenomenon  Solely  a  Hydrodynamic  One? 

Here  we  have  to  do  with  a  body  which,  under  the  effect 
of  the  shock,  is  capable  of  "reacting."  As  a  result  a  chemical 
phenomenon  is  superposed  on  the  hydrodynamic  phenomenon.  Now, 
the  laws  of  chemical  kinetics  are  not  known  sufficiently  well 
under  these  special  conditions  to  permit  us  to  estimate  the 
time  after  which  the  behavior  of  the  explosive  can  no  longer 
be  considered  governed  by  the  laws  of  hydrodynamics  alone. 

We  shall  see  that  In  the  case  of  homogeneous  explosives 
there  exists  an  "Incubation"  time,  a  time  during  which  the 
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chemical  reaction  seems  negligible;  consequently  it  seems 
reasonable  in  this  case  to  treat  the  explosive  as  an  inert 
substance  during  this  period;  it  may  be  considered  that  the 
measurements  carried  out  in  this  way  depend  solely  on  hydro- 
dynamic  phenomena. 

However,  in  the  case  of  heterogeneous  explosives  the 
process  seems  to  be  sharply  different.  The  reaction  is  initi¬ 
ated  almost  immediately  —  which  we  shall  see  later  on,  — 
hence  the  chemical  phenomenon  can  no  longer  be  neglected. 
Nevertheless  the  values  obtained  for  very  wealc  shocks  are 
generally  considered  valid.  Still,  one  is  never  sure  if  the 
shock  has  been  an  inert  one,  and  the  result  which  has  been  ob¬ 
tained  solely  on  the  basis  of  hydrodynamics  should  be  accepted 
with  the  greatest  caution. 

1.4.3.  Shock  Characteristics  of  the  Explosive 

For  reasons  indicated  above,  the  shock  characteristics 
of  solid  granular  explosives  have  been  studied  only  to  a 
limited  extent.  Nevertheless,  for  certain  special  explosives, 
the  methods  used  for  the  determination  of  the  dynamic  adiabatic 
of  solid  inert  substances  —  summarized  above  —  have  been 
used  by  some  researchers  l8]  L 9 3  L 1 2 ]  [16J. 

Another  method,  based  on  the  loss  of  transparence  of 
Plexiglas  (or  glass)  when  subjected  to  a  shock  also  permits, 
after  previous  calibration  under  pressure,  the  measurement  of 
the  shock  velocity  in  an  explosive  sample  of  known  thickness 
placed  between  two  blocks  of  this  substance  [7]. 

To  be  exact,  these  various  methods  do  not  lend  them¬ 
selves  to  a  measurement  of  velocity  but  to  a  measurement  of 
time.  In  order  that  the  average  velocity  Inferred  from  them 
represent  the  physical  phenomenon,  It  is  necessary  that  the 
propagation  of  the  shock  be  constant  during  the  measurement. 

Now  we  know  that  it  is  nothing  of  the  kind,  but  it  Is  impossible 
to  verify  directly  by  these  techniques  If  the  shock  is  reactive 
(acceleration)  or  inert  (decay). 

Hence  a  more  delicate  method  Is  preferable.  V.e  shall 
see  later  on  how  it  Is  possible  to  obtain  a  continuous  diagram 
of  the  course  of  the  shock  In  the  explosive.  Nevertheless  the 
obtainment  of  such  a  diagram  does  not  make  it  possible  to  draw 
any  conclusions  except  in  the  case  where  the  recorded  curve 
is  nearly  linear.  In  the  other  case  we  can  simply  say  that 
the  equations  of  hydrodynamics  do  not  suffice  for  explaining 
the  behavior  of  the  explosive. 
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Chapter  2 

STUDY  OF  THE  INITIATION  OF  DETONATION  BY  SHOCK  Y.AVE 


The  preceding  study  of  shock  generators  has  shown  U9 
that  the  difference  between  the  device  with  barrier  and  the 
Impact  device  consists  essentially  in  the  form  of  the  pressure 
signal  Induced  in  the  receptor.  Hence  we  have  the  right  to 
think  that  the  process  which  leads  to  detonation  is  the  same 
in  both  cases.  Y.hile  a  large  number  of  researchers  have  used 
the  barrier  arrangement,  few  publications  are  available  which 
mention  the  results  obtained  by  impact.  Hence  the  results 
mentioned  below  have,  for  the  most  part,  been  obtained  by  means 
of  experiments  with  barrier. 

2.1.  EXPERIMENTAL  METHODS  USED  FOR  THE  STUDY  OF  THE  GENERATION 
0?  THE  DETONATION  NAVE  IN  THE  EXPLOSIVE 

The  techniques  employed  for  the  determination  of  the 
equations  of  state  of  the  solids  (dynamic  method),  summarized 
above,  have  been  the  first  ones  to  be  employed  by  the 
researchers;  they  have,  nevertheless,  been  soon  dethroned  by 
methods  which  are  better  adapted  to  the  nonstationary  phenomenon 
to  be  studied. 

In  order  to  try  to  determine  the  behavior  of  the  explo¬ 
sive,  the  most  direct  procedure  consists  in  the  continuous 
recording  of  the  progress  of  the  front  of  the  shock  wave  by 
means  of  a  slit  camera. 

Let  us  recall  briefly  the  operation  of  this  type  of 
camera  in  the  version  which  we  have  employed  (rotating  mirror  - 
Brlxner  variant). 

The  image  of  the  phenomenon  to  be  studied  is  formed,  by 
means  of  an  objective,  on  a  slit  placed  in  the  focus  of  a 
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second  objective.  The  parallel  beam  issuing  from  the  latter  is 
reflected  on  a  rotating  mirror.  A  third  objective  forms  the 
image  of  the  slit  on  a  photographic  plate  arranged  in  such  a 
way  that  on  this  plate  space  and  time  correspond  to  two  per¬ 
pendicular  directions. 


Slit  camera  (rotating  mirror). 

(1)  Rotating  mirror;  (2)  Photographic  plate; 

(3)  Slit;  (4)  Object. 

Experimentally,  one  operates  as  follows:  When  the 
initially  transparent  substance  (the  case  of  certain  explosive 
liquids)  becomes  opaque  as  a  result  of  shock,  the  progressive 
darkening  of  a  light  beam  furnished  by  an  argon  flash  placed 
in  the  camera  axis  is  registered  as  a  function  of  time  (record¬ 
ing  by  shadow)  l8,  17,  31 ].  This  method  assumes  that  the 
darkening  immediately  follows  the  shock  front,  which  cannot 
always  be  the  case.  It  nevertheless  makes  it  possible  to 
isolate  the  phenomenon  from  the  edge  effects,  if  the  latter 
do  not  mask  the  phenomenon  along  the  axis  of  the  assembly. 

When  this  procedure  cannot  be  used,  and  if  the  phenomenon 
is  luminous  (the  case  of  a  detonation  wave),  its  evolution  is 
chronologically  recorded  along  a  3ide  parallel  to  the  axis  of 
the  cartridge.  However,  the  one-dimensional  character  is  no 
longer  respected,  and  it  is  no  longer  possible  to  neglect  the 
edge  effects.  In  effect,  the  compatibility  of  pressure  and 
deflection  on  the  faces  of  the  cartridge  with  the  contacting 
medium  Involves  a  lateral  expansion  of  the  medium  which  has 
undergone  shock.  This  perturbs  the  thermodynamic  conditions 
behind  the  shock  front  which  becomes  curved  at  the  vicinity  of 
the  Interface.  Then  the  phenomenon  is  two-dimensional  [7,  12, 
13,  14,  15,  17,  19,  27,  30]. 

To  remedy  this  disadvantage,  another  method  is  usee.  The 
receptor  is  cut  in  the  shape  of  a  wedge  and  the  progress  of  the 
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shock  is  recorded  by  means  of  the  abrupt  variation  which  it 
brings  about  in  the  reflection  of  a  light  beam  of  a  Judiciously 
placed  argon  flash  (the  reflecting  power  of  the  surface  of  the 
wedge  is  generally  obtained  by  means  of  a  metallic  tinsel  a 
few  tens  of  microns  thick). 


f*nt*  (5) 


(1)  Explosive  which  has  been  subjected  to  shock; 

(2)  Argon  flash;  (3)  Fresnel  lens;  (4)  To  slit 
camera;  (5)  Shock  front;  (6)  Inert  explosive. 

The  last  method  presents,  in  addition  to  the  advantage 
of  furnishing  a  continuous  recording  of  the  phenomenon,  that 
of  not  requiring  that  the  latter  be  luminous.  This  is,  no 
doubt,  the  most  precise  procedure  at  the  present  time,  despite 
the  fact  that  one-dimensional  character  of  the  phenomenon  is 
not  strictly  assured  (see  Chapter  3,  paragraph  3*2.1,  and  [9, 
20,  and  29j. 


A  variant  consists  in  placing  a  grid  or  wires  between 
the  flash  and  the  wedge.  The  grid  and  its  image  by  the  tinsel 
are  simultaneously  recorded;  this  image  is  displaced  in  the 
course  of  time  as  a  function  of  the  advance  of  the  shock  front. 
This  method  is  hardly  satisfactory  since,  in  addition  to  the 
fact  that  it  records  the  evolution  of  the  phenomenon  at  a 
finite  number  of  points,  these  points  are  not  fixed  on  the 
receptor. 

Another  method  is  that  of  electric  probes.  It  is  less 
accurate  than  the  preceding  one  since  it  permits  the  recording 
of  the  passage  of  the  shock  only  at  a  finite  number  of  points; 
moreover,  it  requires  an  electronic  setup  with  highly  com¬ 
parable  circuits  so  as  to  prevent  the  introduction  of  disper¬ 
sion  in  the  recording  i6,  9,  17]. 
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(1)  To  08cillo8cope;  (2)  Ketal  probe;  (3)  Schematic  draw¬ 
ing  of  the  principle;  (4)  Schematic  drawing  of  the  assembly 
in  the  receptor  explosive;  (5)  Shock  front;  (6)  Variation 
of  the  front  of  the  electronic  signal. 

Its  operating  principle  is  slightly  different  from  that 
called  upon  when  it  la  used  with  metals.  In  the  latter  case 
the  probe  serves  simply  as  a  circuit  breaker,  where  the  dis¬ 
placement  of  the  metal  closes  a  discharge  circuit.  In  the 
study  of  the  initiation  of  detonation  the  operation  is  based 
on  the  variation  of  electric  resistance  of  the  medium  subjected 
to  shock.  This  variation  is  not  yet  completely  explained.  It 
may  be  said  that  it  depends  roughly  on  the  state  of  advance 
of  the  reaction.  Its  start  is  linked  to  the  passage  of  the 
shock  front,  and  its  variation  in  the  course  of  time  is  a 
function  of  the  conductance  of  the  medium  which  surrounds  the 
two  electrodes  after  the  passage  of  the  shock.  Even  though 
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the  results  obtained  in  this  way  cannot  at  the  present  time  be 
exploited  quantitatively,  they  are  the  only  ones  available 
after  the  oassage  of  the  shock. 

In  view  of  the  impossibility  of  interpreting  the  whole 
recorded  signal,  the  use  of  this  method  reduces,  for  the  moat 
part,  to  the  determination  solely  of  the  moment  of  passage  of 
the  shock  front.  Its  precision  is  quite  low:  in  effect,  since 
the  recorded  signal  depends  on  the  state  of  advance  of  the 
reaction,  its  front  gradually  straightens  out  in  the  course  of 
time  (case  of  heterogeneous  explosives).  Hence  the  recordings 
are  not  directly  comparable  and  for  this  reason  the  exact 
moment  of  the  passage  of  the  shock  front  cannot  be  determined 
with  precision.  On  the  other  hand,  the  introduction  of  elec¬ 
trodes  may  perturb  the  phenomenon  and  lead  to  erroneous  results. 

2.2.  DIFFERENCES  OF  BEHAVIOR  BETWEEN  THE  TWO  TYPES  OF  EXPLOSIVES 
(HOMOGENEOUS  AND  HETEROGENEOUS)  DURING  THE  INITIATION  OF 
DETONATION  BY  SHOCK  WAVE 

The  various  experiments  carried  out  for  the  study  of  the 
initiation  of  detonation  have  revealed  differences  of  behavior 
according  to  the  explosive  under  consideration.  These  differ¬ 
ences  were  found  to  be  considerable,  and  have  led  to  the  clas¬ 
sification  of  the  explosives  into  two  categories: 

--  Homogeneous  explosives  (single  crystals,  liquids); 

--  Heterogeneous  explosives  (granular  solids). 

Let  us  see  how  each  of  these  behaves. 

2.2.1.  Diagram  of  the  Progress  of  the  Shock  Front 

The  study  of  the  detonation  by  shock  wave  of  solid  or 
nonsolid  heterogeneous  explosives  by  one  of  the  above  methods 
L7,  9,  11,  12,  13,  15,  17,  19,  22,  27,  29,  30]  yielded  the 
following  result: 

First  the  shock  front  exhibits  a  progressive  acceleration, 
then  in  a  short  but  measurable  time  it  attains  a  value  close  to 
the  velocity  of  stable  detonation  (see,  however,  Jacobs  and  Seay 
for  some  exceptions  obtained  when  the  density  of  the  charge 
approaches  that  of  the  crystal  l 2 1 ,  28]). 

This  result  is  clearly  different  from  that  obtained  with 
the  homogeneous  explosives  (single  crystals,  liquids).  In 
effect,  in  the  case  of  the  .  atter  it  is  found  that  the  shock 
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first  propagates  at  an  essentially  constant  velocity  up  to  the 
moment  when  the  velocity  suddenly  becomes  greater  than  the 
velocity  of  stable  detonation,  after  which  it  gradually  return® 
to  the  latter  velocity  [8,  20  j. 


•up  I  os  i  f  granutair*  (h4tfrogina)@ 
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(1)  Compression  waves;  (2)  Granular  explosive 
(heterogeneous;;  (3)  Shock;  (4)  Homogeneous 
explosive. 

The  interpretation  which  may  be  given  to  these  two 
recordings  is  as  follows: 

Since  in  the  homogeneous  explosive  the  velocity  of  shock 
propagation  is  constant  (or  slightly  decreasing  --  the  accuracy 
of  measurement  is  low)  during  the  first  part  of  the  phenomenon, 
the  reaction  does  not  Intervene  to  a  notable  extent;  during 
this  period  the  explosive  behaves  as  an  inert  substance.  Then, 
abruptly,  the  reaction  is  initiated  at  the  interface  —  the 
region  which  was  the  first  to  undergo  the  shock.  This  reaction 
Immediately  leads  to  detonation.  The  latter  is  then  propagated 
In  the  compressed  explosive.  Its  velocity  is  greater  than  the 
normal  detonation  velocity  of  the  explnolve.  It  overtakes  the 
shock  wave  and,  for  some  time,  imposes  a  superdetonatlon  on 
the  nonshocked  explosive.  The  final  evolution  leads  to  the 
steady  state  of  detonation  [8,  10]. 

By  contrast,  In  the  case  of  the  heterogeneous  explosive, 
the  reaction  is  triggered  Immediately  and  it  progressively 
accelerates  the  shock  front.  It  then  follows  that  at  no  time 
can  the  explosive  be  considered  as  inert.  At  a  given  moment  a 
rapid  velocity  Jump  leadB  to  the  steady  state  which  is  attained 
asymptotically  [93. 
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It  la  nevertheless  possible,  for  each  of  the  two  types 
of  explosive,  to  determine  a  time  which  we  shall  call  "Induc¬ 
tion  period"  after  which  the  steady  state  of  detonation  Is 
established.  This  parameter  will  be  rather  poorly  defined  In 
the  two  cases,  since  It  Is  found  that  It  Is  not  obtained  after 
an  abrupt  discontinuity  but  attained  gradually  by  a  decelera¬ 
tion  of  the  shock  front  In  the  homogeneous  case,  and  by  an 
acceleration  of  the  shock  front  In  the  heterogeneous  case. 

Even  tnough  Its  absolute  measurement  Is  not  accurate, 

It  Is  nevertheless  possible  to  study  Its  relative  variation  as 
a  function  of  the  Initial  conditions. 

2.2.2.  Threshold  of  Initiation  of  the  Detonation 

If  the  shock  Is  too  weak  the  signal  obviously  decays 
without  causing  detonation.  Hence  It  is  only  above  a  certain 
threshold  value  that  the  chemical  process  which  leads  to 
detonation  Is  initiated.  This  threshold  Is  very  different 
according  to  the  explosive  being  tested.  Whereas  a  few  kilo- 
bars  suffice  to  detonate  the  heterogeneous  explosives,  a  large 
signal  is  necessary  for  bringing  about  the  detonation  of  homo¬ 
geneous  signals  (85  kb  for  nltromethane ,  112  kb  for  a  pentrlto 
crystal)  l20,  29]. 

2.2.3.  Initial  Point  of  the  Start  of  Detonation 

In  the  case  of  granular  heterogeneous  explosives,  the 
experiments  have  been  carried  out  on  cartridges  of  rather 
small  dimensions.  In  this  case  the  "edge  effects"  intervene 
and  the  phenomenon  no  longer  exhibits  the  one-dlnen3ional 
characteristic . 

These  experiments  have  nevertheless  permitted  to  demon¬ 
strate  the  fact  that  the  detonation  does  not  take  place 

--  In  time:  immediately  after  the  shock; 

--  In  space:  at  the  Interface. 

It  should  be  mentioned  that  this  result  does  not  contra¬ 
dict  that  mentioned  above  for  the  same  explosives.  In  effect, 
while  the  chemical  reactions  Indeed  take  place  at  the  Inter¬ 
face,  the  shock  becomes  gradually  "more  and  more  reactive"  but 
nevertheless  it  cannot  be  said  tnat  detonation  has  occurred. 

These  results  were  first  published  by  Hertzberg  and 
Walker  l 1 9 ]  --  initiation  of  detonation  by  jeans  of  a 
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detonator  --  and  have  since  been  verified  by  numerous  authors, 
using  different  methods: 

—  The  Integral-image  camera  (10^  images  per  second) 
clearly  shows  that  the  Initially  luminous  point  in  the  receptor 
1 8  situated  at  a  certain  distance  from  the  interface; 

—  The  slit  camera,  during  a  chronological  study  of  the 
luminosity  of  the  edge  of  the  cartridge,  registers  this  lag  in 
time  and  space  in  a  quantitative  manner; 

--  The  method  of  probes  reveals  a  progressive  conduc¬ 
tivity  of  the  medium  after  the  passage  of  the  shock,  In  direct 
proportion  to  the  distance  from  the  interface  [9,  17,  2?]. 

Beginning  at  this  first  point  of  initiation,  a  detonation 
wave  is  propagated  in  the  inert  explosive.  Another  detonation 
wave  Issuing  from  the  same  point  may  rise  toward  the  Interface 
if  the  already  shociced  explosive  has  reacted  only  to  a  slight 
extent.  This  phenomenon,  called  "retonation , "  has  been  ob¬ 
served  experimentally  1.11,  12,  1 3 »  17,  27  J .  Its  interpretation 
has  generally  been  that  it  is  a  consequence  of  the  expansion 
effects  due  to  the  edges  of  the  cartridge,  since  the  one- 
dimenslonpl  character  of  the  phenomenon  was  no  longer  assured. 

By  contrast,  in  the  case  of  homogeneous  explosives,  the 
photos  taken  by  means  of  the  integral- image  camera  [8]  show 
that  a  clearly  detached  shock  precedes  the  detonation.  The 
slit  camera  [8,  31 ]  permits  the  recording  of  a  rather  low 
luminosity  --  associated  with  the  detonation  wave  in  the  com¬ 
pressed  explosive  --  prior  to  the  much  more  Intense  luminosity 
of  the  state  of  superdetonatlon  in  the  nonshockea  explosive. 
Finally,  measurements  by  electric  probes  [8,  31 J  show  that 
after  a  certain  "incubation"  period  the  reaction  is  abruptly 
initiated  at  the  interface. 

These  various  experiments  do  not  furnish  a  direct  meas¬ 
urement  of  the  first  point  of  the  initiation  of  the  detonation, 
but  they  constitute  excellent  indirect  proofs  with  which  to 
confirm  the  idea  that  the  detonation  is  initiated 

—  In  the  case  of  heterogeneous  explosives:  within  the 
explosive; 


--  In  the  case  of  homogeneous  explosives:  at  the  inter¬ 
face. 
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2.2.4.  Influence  of  the  Intensity  of  the  Pressure  Signal 

In  proportion  to  the  increase  of  the  Intensity  of  the 
pressure  signal,  the  steady  3tate  of  detonation  sets  in  sooner 
in  time  and  closer  to  the  Interface  in  space.  Kowpver  these 
variations  are  notably  different  for  the  two  types  of  explosive. 

V/hlle  in  the  case  of  the  heterogeneous  explosives  there 
is  a  progressive  variation  as  a  function  of  the  intensity  of 
this  signal,  this  Is  not  the  case  for  the  homogeneous  explosives 
where,  once  the  threshold  has  been  reached,  the  variation  is 
abrupt.  In  the  case  of  nltromethane,  for  example  l8]>  when  the 
shock  changes  from  86  kb  to  89  kb  (a  3. 3/2  increase),  the  "Induc¬ 
tion"  time  is  decreased  by  26,'i. 
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(l)  Distance  in  mm;  (2)  Pressure  in  kb; 

(3)  Extracted  from  l 9 3 . 

This  comparison  leads  one  to  think  that  In 
the  case  of  the  homogeneous  explosives  there  exists  a  "threshold 
state"  below  which  the  process  leading  to  detonation  is  not 
initiated.  However  once  this  threshold  has  been  exceeded  the 
phenomenon  is  very  abrupt.  In  the  case  of  the  heterogeneous 
explosives,  on  the  other  hand,  the  duration  of  the  transitory 
state  varies  continuously,  and  assures  a  more  gradual  construc¬ 
tion  of  the  detonation  wave. 

2.2.5.  Effect  of  the  Initial  Temperature 

The  initial  temperature  also  Intervenes  in  a  clear 
manner  in  the  differentiation  of  these  two  tjpes  of  explosive. 

In  the  case  of  the  homogeneous  explosive,  a  temperature  in¬ 
crease  of  30°K  divides  by  3  the  time  necessary  for  the  estab¬ 
lishment  of  the  steady  state  of  detonation  [8j  while  in  the 
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case  of  the  heterogeneous  explosive,  the  influence  is  much 
less,  and  undetectable  for  such  a  small  temperature  variation 
[9]. 


2.3.  MECHANISMS  OF  THE  INITIATION  OF  DETONATION 

Various  mechanisms  have  been  proposed  to  explain  the 
physical-chemical  phenomena  which  take  place  prior  to  the 
establishment  of  the  steady  state  of  detonation.  Few  of 
these  mechanisms  hove  received  any  experimental  support  in 
view  of  the  microscopic  scale  on  which  these  mechanisms  take 
place,  and  also  because  of  the  fact  that  at  the  present  time 
it  is  only  possible  to  study  experimentally  the  macroscopic 
aspect  of  this  problem. 

2.3*1 •  The  Hot  Spots 

A  frequently  proposed  mechanism  (cf.  Bowden  and  co¬ 
workers)  is  the  creation,  within  the  explosive,  of  hot  spots 
by  adiabatic  compression  of  the  occluded  gaseB. 

Although  this  hypothesis  has  received  satisfactory 
experimental  support  in  the  case  of  shocks  cf  low  intensity, 
this  is  not  so  in  the  case  which  is  of  interest  to  us.  It 
permits  the  qualitative  explanation  of  certain  well-known 
facts.  For  example  we  know  that  the  lower  the  charge  density 
of  a  granular  explosive,  the  easier  it  is  to  initiate  the 
detonation.  Likewise  the  fine-grain  charges  are  more  sensitive 
to  shock,  and  this  difference  of  behavior  as  a  function  of 
grain  size  gradually  becomes  attenuated  as  one  approaches  the 
density  of  the  crystal.  Hence  the  predominant  role  seems  to 
be  played  by  the  occluded  gases. 

If  the  phenomenon  is  a  thermal  one,  then  the  kind  of 
gas  1 8  an  Important  factor.  In  effect,  depending  on  the  value 
of  the  ratio  of  specific  heats  of  this  gas,  the  temperature 
at  the  end  of  the  compression  may  vary  considerably  for  a 
shock  of  the  same  intensity. 

Experiments  carried  out  with  argon  (7  =  1.67),  methane 

(7  =  I.31)  —  the  case  of  a  perfect  gas  leads  to  t,  •  t, 

and  l_i_l  -  0.40  for  7  =  1.67  while  it  is  0.246  when  7  =  1 . 3 1  — 
y 

have  shown  that  while  the  discounted  variation  indeed  takes 
place  in  the  anticipated  direction,  it  is  very  far  from  the 
calculated  value.  Moreover,  by  replacing  the  gas  bubbles  with 
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balls  made  of  different  materials  —  particularly  tungsten  — 
the  obtained  delays  are  of  the  same  order  of  magnitude.  The 
situation  is  the  same  when  a  hieh  vacuum  is  created  [9,  29 J. 

Hence  these  various  experiments  make  us  believe  that 
the  Influence  of  the  hot  spots  cannot  be  considered  as  a  funda¬ 
mental  mechanism,  since  the  Initiation  of  the  detonation  prac¬ 
tically  does  not  involve  the  temperature  after  the  passage  of 
the  shock  wave. 

Nevertheless,  both  in  the  case  of  the  homogeneous 
explosives  and  that  of  heterogeneous  explosives  it  is  well 
observed  that  the  detonation  starts  out  from  these  points. 

This  cause  and  effect  relationship  has  been  demonstrated 
in  several  ways: 

In  the  case  of  the  homogeneous  explosives,  the  intro¬ 
duction  of  fine  gas  bubbles  at  the  interface,  or  the  scratching 
of  grooves  on  the  barrier  [8,  31 ]»  has  shown  that  the  reaction 
starts  from  these  bubbles  or  grooves,  provided  that,  in  the 
case  of  the  bubbles,  their  diameter  is  not  less  than  a  certain 
minimum  value  (approx.  0.7  mm  for  nitromethane ) . 

In  the  case  of  the  heterogeneous  explosives,  realization 
of  these  fine  grooves  on  the  surface  of  the  explosive  subjected 
to  the  shock  wave  considerably  decreases  the  delay  which  pre¬ 
cedes  the  establishment  of  the  stabJe  detonation  regime  [9; 

31,  P  56]. 

Mader  has  proposed  that  the  experimental  results  obtained 
with  homogeneous  explosives  be  interpreted  by  taking  into  ac¬ 
count  the  hydrodynamic  character  of  the  phenomenon  ("hydro- 
dynamic  hot  spot")  l 26 3  • 

Basing  himself  on  the  theoretical  results  of  Evans, 

Harlow  and  Kelxner  L 1 6 ]  who  have  calculated  that  the  inter¬ 
action  of  a  shock  with  a  bubble  created  within  the  perturbed 
medium  a  hot  spot  which  possesses  essentially  the  same  volume 
as  the  bubble,  he  confirmed  theoretically  the  dimensional 
effect  obtained  experimentally. 

The  "hydrodynamic  hot  spot"  has  a  temperature  (case  of 
the  "temperature  hot  spot")  and  a  pressure  (case  of  the  "pres¬ 
sure  hot  spot")  which  are  above  those  of  the  environment. 

The  law  of  chemical  kinetics  employed  is  Arrhenius'  law.  It 
can  be  seen  then  that,  depending  on  the  diameter  of  the  sphere 
constituting  the  hot  spot,  the  reaction  may  or  may  not  have 
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the  time  to  be  Initiated  before  the  expansion  waves  propagating 
toward  the  center  of  the  hot  point  had  sufficiently  decreased 
the  latter's  temperature. 

The  transposition  of  this  study  to  the  heterogeneous 
medium  can  only  be  suggested  at  the  present  time.  But  It  is 
hard  to  see  how  we  will  be  able  to  explain,  by  means  of  a  model 
which  Imposes  a  minimum  dimension  of  the  hot  spot,  why  the  fine¬ 
grained  explosives  have  the  highest  sensitivity.  The  minimum 
intensity  of  the  shock  necessary  to  Initiate  the  detonation 
Is  much  lower  (a  few  kllobars  Instead  of  80-100  kb),  hence  the 
temperature  of  the  hot  spot  is  much  lower,  and  hence  Arrhenius' 
law  calls  for  more  time  to  bring  about  the  initiation,  whereas 
it  Is  found  experimentally  that  there  is  no  delay  due  to  "incu¬ 
bation"  in  the  case  of  these  explosives,  in  contrast  with  the 
homogeneous  explosives. 

On  the  other  hand,  since  Arrhenius'  law  entails  an  abrupt 
release  of  detonation,  the  gradual  acceleration  of  the  reactive 
shock  front  of  the  granular  explosive  seems  difficult  to  explain. 

2.3.2.  Other  Mechanisms 

By  contrast,  the  above-mentioned  experiments  support  the 
fact  that  If  a  discontinuity  exists  within  the  explosive,  the 
detonation  19  initiated  from  that  point.  This  discontinuity 
may  be  present  In  various  forms:  gas  bubbles,  foreign  bodies, 
grooves,  etc.  The  physical  fact  which  may  be  associated  with 
It  Is  the  perturbation  which  It  entails  with  respect  to  the 
shock  front.  The  latter  Is  no  longer  plane.  Various  mechanisms 
have  been  proposed  on  the  basis  of  this  remark:  the  formation 
of  these  sinusoid  shock  waves  may  lead  to  the  pulverization  of 
the  grains  of  explosive.  This  is  so  since,  if  they  are  con¬ 
vergent,  they  produce  local  overpressures  and  bring  about 
hollow-charge  effects;  if  they  are  divergent,  they  produce  In 
the  grains  tensions  which  may  lead  to  mechanical  ruptures 
favored  by  the  crystal  defects  (scaling)  or  chemical  ruptures 
due  to  the  changes  In  the  molecular  bonds  l3*  4,  5»  7,  28].  An 
Increase  in  temperature  may  also  produce  bursting  of  the  grain 
as  a  result  of  expansion. 

All  this  leads  —  as  suggested  by  Andreev  [2]  --  to  the 
formation  of  a  fine  suspension  whose  pressure  is  abruptly  In¬ 
creased  by  the  explosion,  and  this  pressure  increase  assures 
the  gradual  generation  of  the  stable  detonation  regime,  how¬ 
ever,  the  assumption  that  this  mist  acts  by  impact  or  by  fric¬ 
tion  on  the  following  grains  does  not  represent  a  supplementary 
explanation  but  takes  up  on  the  microscopic  scale  the  mechanisms 
which  have  been  refuted  ma cro scopi cally . 
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2.4 


SPECIFIC  RESULTS  OBTAINED  BY  IMPACT 


i 


The  only  known  results  are  those  of  Brown  and  Whitbread 
[6],  Although  they  are  not  one-dimensional,  they  nevertheless 
permit  the  estimation  of  the  influence  of  the  special  parameters 
linked  to  this  mode  of  generation  of  the  shock  wave.  The  pres¬ 
sure  signal  obtained  by  this  procedure  is,  in  effect,  of  a 
rectangular  shape,  constant  during  a  certain  period  of  time  and 
then  decreases  rapidly  (see  Chapter  on  the  study  of  shock  gen¬ 
erators).  This  particular  form  of  the  pressure  signal  has  made 
it  possible  to  show  that,  in  order  to  bring  about  the  detona¬ 
tion,  its  intensity  and  its  time  of  application  must  be  greater 
than  the  minimum  values  pm  and  Tm  (t^  being  the  minimum  associ¬ 
ated  with  pm). 

If  the  Intensity  is  less  than  this  minimum  (Pm)>  the 

signal,  regardless  of  its  duration,  cannot  assure  the  detona¬ 
tion  of  the  explosive.  On  the  otr.er  hand,  if  the  intensity  is 
greater  than  this  minimum,  it  is  possible  that  the  detonation 
will  be  established,  even  if  the  time  of  application  is  less 
than  the  minimum  time  Tm. 

2.5.  ANOTHER  MODE  OF  GENERATION  OF  DETONATION:  THE  TRANSITION: 

DEFLAGRATION  -  DETONATION 

The  transition:  deflagration  -  detonation,  which  to  be 
sure  is  slower  than  the  transition:  shock  -  detonation,  is 
nevertheless  capable  of  furnishing  certain  data  with  regard  to 
the  mechanism  leading  to  the  initiation  of  detonation  in  the 
explosive. 

Whereas  the  steady  states  of  deflagration  and  detonation 
are  sufficiently  well  known,  the  passage  from  one  to  the  other 
is  much  less  well  known,  desnlte  numerous  studies  il8»  23,  24, 
25,  32,  33 J . 

It  is  nevertheless  obvious  that  the  establishment  of  the 
steady  detonation  state  by  this  process  is  necessarilv  accom¬ 
panied  by  the  formation  of  a  shock.  Various  authors  l23,  25, 

33 J  believe  that  the  formation  of  this  shock  is  the  direct 
cause  leading  to  the  detonation.  By  means  of  this  hypothesis, 
the  transition  may  be  described  as  follows:  The  rapid  Increase 
of  pressure  behind  the  combustion  front  produces  compression 
waves  which  are  propagated  in  the  non-burned  explosive  in 
front  of  the  flame.  The  combining  of  these  elementary  waves 
produces  a  shock  wave  which  initiates  the  detonation. 


32 


Experimental  measurements  made  by  Macek  [24]  show  well 
the  exponential  variation  of  the  pressure  at  the  shock  front 
with  time  during  the  transition:  deflagration  -  detonation,  and 
the  theoretical  study  carried  out  by  Zovko  and  Macek  L 33 J  with 
this  hypothesis  leads  to  a  rather  satisfactory  representation 
of  the  whole  phenomenon. 

The  various  studies  which  have  been  carried  out  tend  to 
prove  that  the  pressure  is  indeed  the  predominant  factor  in  the 
establishment  of  the  detonation.  Thus  the  assumption  that  it 
is  also  the  predominant  factor  in  defining  the  reaction  rate  of 
the  explosive  is  not  without  foundation. 
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Chapter  3 

STUDY  OF  THE  EXPERIMENTAL  APPARATUS  EMPLOYED 


3.1.  THE  SHOCK  GENERATOR 

3.1.1.  The  Lifting  of  a  Plate 

The  experimental  procedure  chosen  for  bringing  about 
the  Impact  Is  that  of  lateral  projection,  by  the  so-called 
"plate  lifting." 

As  we  shall  recall,  it  consists  in  projecting  a  metal 
lining  of  low  thickness  by  means  of  an  explosive  initiated  in 
such  a  manner  that  its  detonation  front  is  perpendicular  to 
the  plate  Experience  shows  that  after  a  zone  of  accleration, 
the  plate  becomes  essentially  plane  and  assumes  a  direction 
which  makes  a  constant  angle  with  that  which  it  had  originally 
[9]  (Plate  I,  Fig.  1). 

The  shock-  and  wave- reflection  phenomena  induced  in  the 
lining  by  the  detonation  of  the  explosive  are  sufficiently 
fast  to  be  able  to  neglect  the  zone  of  acceleration  in  first 
approximation,  and  assume  that  the  plate  has  been  set  into 
motion  instantaneously. 

nils  may  be  represented  schematically  as  follows: 

During  a  unit  tine  Interval  during  which  the  detonation  wave 
has  shifted  from  point  M  to  point  M'  (Ml!1  =  D,  detonation 
velocity  of  the  explosive),  the  particles  of  the  material 
originally  situated  at  M  have  reached  P. 

•4 

The  velocity  vector  V  of  the  particle  is  defined  by 
ftf.  1)  *  |  -  |  V  .  3  D  .in  | 
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and  the  plate  hae  an  overall  movement  defined  by  the  vector  U 

(6,  tf)  •  |  -  9  0  •  D  *in  f 


Hence  the  plat^e  haa  an  overall  displacement  in  the 
direction  of  vector  if  (component  of  vector  V  normal  to  the  plane 
of  the  raised  plate  )_^by  gliding  along  the  tangential  component 

of  this  same  vector  V. 

To  be  sure,  the  friction  phenomena  studied  by  Bowden 
[2,  3]  at  lesser  velocities  nevertheless  show  that  the  delays 
necessary  for  assuring  the  detonation  of  the  explosive  are  very 
much  greater  than  those  which  we  wish  to  study  here.  Hence  it 
seems  that  these  delays,  while  not  negligible  in  absolute 
value  --  the  velocity  of  lateral  displacement  is  equal  to 
2  D  sin2  <p/2  —  do  not  have  the  time  to  perturb  the  principal 
phenomenon  linked  to  the  overall  velocity  U  of  the  plate,  and 
we  shall  assume  that  it  is  this  value  that  represents  the 
velocity  of  Impact  of  the  plate  projected  on  the  target. 

3.1.2.  Practical  Execution 

The  experimental  apparatus  comprises  three  parts: 

—  The  primer  which  assures  the  detonation  of  the 
explosive ; 


—  The  explosive  in  the  form  of  a  thin  rectangular 

plate; 


--  The  metallic  plate  to  be  projected,  in  contact  with 
the  explosive. 
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To  obtain  a  plane  two-dimensional  phenomenon,  it  is 
necessary  to  assure  the  simultaneous  detonation  of  one  of  the 
ends  of  the  explosive  plate.  Since  the  latter  is  thin  com¬ 
pared  to  its  other  two  dimensions,  the  problem  becomes  that  of 
bringing  about  the  simultaneous  detonation  along  a  straight- 
line  portion.  The  primer  Is  then  constituted  of  a  "linear 
generator."  The  model  retained,  as  described  by  Erkman  [6], 
is  a  surface  constituted  of  a  portion  of  a  cone  of  revolution 
and  a  plane. 


Schematic  drawing  of  a  linear  generator. 

(l)  -  Lengths  A,  are  all  equal  to 

AB . 

Such  "surface^"  ha^e  been  molded  from  explosives 
(average  thickness  10  mm)  and  the  simultaneousness  of  the 
arrival  of  the  detonation  wave  on  the  exit  surface  was  checked 
by  means  oi  a  silt  camera.  The  maximum  error  Is  equal  to 
♦  0.07  psec. 

The  explosive  plate  Is  sufficiently  wide  so  that  the 
edges  of  the  lining,  which  are  raised  to  a  lesser  extent  than 
the  central  part  since  they  do  not  undergo  the  same  thrust  as 
a  result  of  the  lateral  expansions  to  which  the  detonation 
products  are  subjected,  do  not  perturb  the  two-dimensional 
character  of  the  phenomenon.  They  are  sufficiently  long  so 
that  the  states  of  stable  detonation  and  constant  life  are 
established . 

Through  these  precautions  a  "useful"  zone,  measuring 
several  centimeters  on  a  side,  could  be  obtained.  Tests  for 
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planeness  and  for  simultaneity  on  Impact  --  have  been  carried 
out  In  the  two  dimensions  by  means  of  the  silt  camera,  fhe 
useful  part  Is  plane  to  within  4  0.05  fjsec,  or  for  an  Impact 
velocity  of  1,000  m/sec,  to  within  4  50  p.  This  Is  the  mech¬ 
anical  precision  with  which  the  thickness  of  the  lining  Is 
defined. 

3.1.3*  Experimental  Measurement  of  the  Impact  Velocity  of  the 

Projected  Plate 

The  determination  of  U  =  D  sin  cp  depends  on  the  exper- 
mental  measurement  of  the  two  quantities: 

--  D,  the  velocity  of  detonation  of  the  explosive; 

--  f,  the  angle  of  lift  of  the  projected  lining. 

The  velocity  of  detonation  of  the  explosive  was  obtained 
by  the  classic  method  consisting  in  measurinc  the  time  elapsed 
during  the  passage  of  the  detonation  wave  between  two  points 
materialized  by  probes.  Their  positioning  (4  0.5  mm  for  a  base 
of  100  mm)  combined  with  the  use  of  an  electronic  chronometer 
having  an  accuracy  of  4  1/40  pisec  leads  to  a  measurement  of  the 
velocity  of  detonation  having  an  accuracy  of  the  order  of  1,». 

For  the  measurement  of  angle  >f',  Its  order  of  magnitude 
may  be  obtained  in  various  ways: 

--  Photographically  by  means  of  the  Integral-Image 
camera  (10^  lmages/sec); 

—  By  flash  radiography  (exposure  time  aoout  0.1  ^sec). 

As  a  result  of  the  inaccuracy  regarding  the  exact 
position  of  the  lining  on  these  recordings  (blur  due  to  the 
non-negllglble  exposure  time  and  to  the  main  of  the  film),  the 
measurement1*  carried  out  do  not  permit  an  accuracy  meater 
than  4  30'  In  the  measurement  of  angle  ^  . 

In  order  to  obtain  a  better  result,  the  slit  camera 
was  used.  The  schematic  drawing  of  the  assembly  Is  shown 
below:  Lnext  pagej 

Angle  <t,  whose  magnitude  is  close  to  angle  f,  is  ad¬ 
justed  prior  to  the  experiment  by  means  of  an  autocollimator 
device  for  angle  measurement.  A  series  of  standard  wedges 
make  It  posslole  to  obtain  a  precision  of  the  order  of  one 
minute.  However,  such  a  precision  is  an  illusory  one,  since 
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despite  the  fact  that  the  shock-generating  assembly  is  rigid, 
the  metal  lining  to  be  projected  does  not  have  a  planeness 
which  is  sufficiently  stable  mechanically  to  guarantee  this 
precision.  However  an  accuracy  to  within  a  few  minutes  mev  be 
considered  reasonable. 


Schematic  drawing  of  the  assembly  for  the 
measurement  of  *f. 

(1)  -  plate  of  explosive;  (2)  -  lining; 

(3)  -  target  device;  (4)  direction  of  ob¬ 
servation  . 

The  projected  plate  strikes  a  target  which  lights  up 
under  the  Impact.  This  luminosity  is  obtained  by  an  argon 
chamber.  The  slit  camera  registers  this  phenomenon  chrono¬ 
logically  . 

The  treatment  of  these  firings  Is  carried  out  in  the 
following  manner:  The  preliminary  recording  on  film  of  two 
luminous  points  placed  in  the  exact  area  where  the  target  will 
be  placed  permits  the  measurement  of  the  magnification  of  the 
whole  setup  (k  =  Xj/x),  where  x1  is  the  distance  between  the 

two  lines  recorded  on  the  film,  and  x  the  distance  between  the 
two  luminous  points. 

The  measurement  of  tne  slope  (tan  V  )  of  the  signal 
recorded  on  the  film  ( V'  is  the  angle  between  this  trace  and 
the  line  perpendicular  to  the  time  axis)  a9  well  as  the  sweep¬ 
ing  velocity  (V)  of  the  camera  at  the  moment  of  firing  then 
permit  the  determination  of  <?’  by  the  equation 


tg  *> 


cos  a  +  k  — tg  T 
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,T 
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Tan  r  la  positive  if  cp  is  smaller  than  <x.  ,  and  negative  in 
the  opposite  case. 


Adjustment  of  angle  cc  . 

(1)  -  autocollimator;  (2)  -  standard 
wedge;  (3)  -  lining;  (4)  -  target. 

The  accuracy  of  the  measurement  of  P  is  not  excellent, 
since  the  signal  rect rded  is  not  strictly  linear  (deviation 
♦  0.0b  ^isec).  Nevertheless  it  seems  that  the  maximum  error  is 
about  1/2  degree,  which  causes  an  error  in  the  value  of  <p  of 
less  than  5  minutes.  Finally,  the  angle  of  lift  of  the  plate 
is  determined  to  within  ♦  10'. 

3.1. 4.  Theoretical  Aspect  of  the  Lateral  Projection 

Richter  L6],  has  studied  in  1945.  by  classical  mechanics, 
the  movement  of  the  metal  lining  subjected  to  the  effect  of  the 
explosive.  He  made  certain  simplifying  assumptions  regarding 
the  behavior  of  the  products  of  detonation  and  of  the  lining, 
calling  upon  fluid  mechanics  (theory  of  shock  waves  and  method 
of  characteristics)  only  for  a  qualitative  explanation  of  the 
phenomena.  These  assumptions  enabled  him  to  give  this  problem 
simple  and  completely  integrable  solution. 

Moreover,  he  made  an  effort  to  determine  the  effect  of 
various  factors  on  the  idealized  solution.  This  has  enabled 
him  to  show  that  the  compressibility  of  the  metal,  its  resistance 
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to  deformation  and  the  pressure  gradients  appearing  in  its 
interior  only  change  the  idealized  solution  by  a  few  percent. 

After  a  brief  recapitulation  of  Richter's  theory  adapted 
to  the  two-dimensional  case,  we  shall  treat  the  same  problem 
by  the  method  of  characteristics. 

—  Richter's  Method: 

The  metal  plate  is  considered  as  being  constituted  of  a 
sequence  of  independent,  Juxtaposed  material  points.  Since  the 
flow  is  steady,  we  look  for  the  path  of  one  of  these  points 
with  which  mass  £,  is  associated,  which  is  the  eurface  density 

of  the  lining  where  is  the  density  of  the  metal  employed  in 
the  form  of  a  plate  of  thickness  fc. 

In  the  reference  system  linked  with  the  detonation  front, 
the  acceleration  of  such  a  point  M  is  given  by  the  two  equations 


Y»  *  0  tangential  acceleration  (constant 
v,  velocity) 

t.  ■-J-  normal  acceleration 


Hence  this  point  hau  a  movement  defined  by  the  sole 
equation  p.  e  y,  *  p,  *  p,  where  p?  represents,  at  this  point,  the 

pressure  of  the  detonation  products,  and  P2  the  pressure  of  the 
medium  adjacent  to  the  plate. 


(1)  -  detonation  products;  (2)  -  solid 
explosive. 
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Since  the  flow  Is  assumed  to  be  steady,  v  Is  equal  to 
D  (=  ds/dt),  s  Is  the  curvilinear  abscissa  of  point  K,  measured 
for  example  from  the  detonation  front). 

If  ©  is  the  angle  between  the  tangent  to  the  path  of 
point  M  and  the  reference  axis  (separating  the  explosive  and 
the  metal  prior  to  detonation),  R  =  ds/d©  and  the  preceding 
equation  is  transformed  into 


p.  1 


Pi 


p. 


(D 


The  assumptions  made  by  Richter  with  regard  to  p^  and 
P2  are  as  follows: 

To  determine  p i ( s  ) ,  he  associates  two  expansions:  the 
first  produced  by  the  lifting  of  the  plate,  the  second  due  to 
the  finite  thickness  of  the  explosive. 

If  the  thickness  of  the  explosive  is  infinite,  only  the 
first  expansion  exists.  The  latter  may  be  estimated  by  means 
of  the  p(0)  curve  derived  from  busemann's  epicycloid  by  assuming 
the  existence  of  an  lsentropic  expansion  at  P  which  is  con¬ 
stant  and  equal  to  3.  Richter  likens  this  curve  to  a  straight 
line  whose  initial  conditions  are  as  follows: 

9  ■  <p. - »  p,  *  o 

9-o - -  Pj  .  p  (detonation  pressure) 


(1)  -  expansion  curve  of  the  detonation 
products;  (2)  -  "heart-shaped"  curve  of 
adjacent  medium;  (3)  -  Richter's  straight 

line. 
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The  angle  of  the  smokes  <fm  Is  the  angle  of  deviation  of 

the  detonation  products  when  the  explosive  detonates  without 
lining  (£=  0).  V.e  obtain, 


whence 


*  ('  -T.) 


1  £2 
%  «*■ 


When  the  thickness  of  the  explosive  is  finite  (=  e), 
there  is  added  to  this  expansion  the  wave  bundle  centered  on 
the  edge  of  the  detonation  front,  on  the  side  opposite  to  the 
plate.  Richter  evaluates  it  for  this  two-dimensional  case  as 

.  .  A  Pi  where  /'■is  a  parameter  which  depends  solely  on  the 
explosive. 


The  combination  of  these  two  expansions  leads  to 


p  de 


(2) 


The  pressure  p2  is  furnished  by  the  determination  of  the 
supersonic  flow  of  the  medium  adjacent  to  the  lining  (in  gen¬ 
eral,  air).  It  Is  of  the  order  of  a  few  hundred  kg/ cm2  only, 
and  cannot  modify  the  behavior  of  the  lining.  Hence  it  is 
legitimate  to  neglect  it. 

Combination  of  equations  (1)  and  (2)  leads,  after  inte¬ 
gration,  to 


9 


e-c») 


where  C  is  given  by 

c _ e _ .4  -L 

p.  t  D'  De 


When  s  approaches  infinity,  ©  approaches  <p,  the  angle 
of  lift  of  the  metal  plate. 

i  1  DX.  p  c 

v’T*~r  -t- 
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The  experimental  result9  which  we  have  obtained  by  the 
method  presented  In  the  preceding  paragraph  confirm  this  linear 

law  ^  (p.  in  the  range  of  values  of  £  and  e  employed 

(Plate  I,  Fig.  2). 

-  Method  of  Characteristics 

In  order  to  Improve  the  theoretical  determination  of 
the  lift  of  the  plate,  it  is  necessary  to  follow  in  a  more 
precise  manner  the  behavior  of  the  different  media:  detonation 
products,  lining. 

In  effect,  Richter* s  theory  leads  to  satisfactory 
results  if  the  respective  values  of  £  and  e  are  sufficiently 
small  (a  few  millimeters  for  £,  a  few  centimeters  for  e). 
Nevertheless,  it  furnishes  a  dimensionless  result  in  e/e,  but 
it  1 8  obvious  that  there  exists  a  minimum  value  of  e  beyond 
which  the  explosive  no  longer  intervenes  to  raise  the  metal 
plate.  Likewise  the  value  of  the  lifting  angle  obtained  for  a 
large  £  is  not  confirmed  experimentally.  This  has  to  do  no 
doubt  with  the  assumptions  made  with  regard  to  both  the  laws 
of  expansion  of  the  detonation  products  and  the  fact  that  it 
is  only  the  density  which  intervenes  in  the  characterization 
of  the  material  which  constitutes  the  lining. 

The  theoretical  study  by  means  of  the  method  of  charac¬ 
teristics  Itself  requires  that  certain  assumptions  be  made. 

The  most  important  of  these  assumptions  —  since  we  use 
the  general  hydrodynamic  equations  without  a  conduction  or 
viscosity  term  --  is  that  the  metal  is  a  fluid  which  behaves 
as  such  at  the  pressures  under  consideration.  This  is  quite 
Justified,  as  is  shown  by  current  studies  of  its  behavior  at 
high  pressure  Li,  7 J .  The  plastic  character  is  no  longer 
taken  into  consideration,  but  an  analysis  of  its  influence 
LP.  Beatrix,  unpublished  results]  has  revealed  that  this 
characteristic  modified  the  behavior  of  the  metal  only  to  a 
negligible  extent. 

Experimentally  the  explosive  and  the  metals  are  sur¬ 
rounded,  except  under  special  conditions,  by  air  at  atmospheric 
pressure,  be  have  already  seen  that  the  influence  of  the 
latter  is  negligible.  Theoretically  we  have  not  taken  it  into 
consideration.  Hence  this  amounts  to  considering  that  physi¬ 
cally  the  experiment  is  carried  out  in  vacuo. 


It  is  not  our  intention  to  present  at  this  point  the 
mathematical  theory  of  characteristics  (cf.  [4]  and  L 5 3  in  this 
connection),  rather,  we  shall  be  content  with  pointing  out  the 
various  particular  aspects  of  its  use  in  the  problem  of  the 
lifting  of  the  plate. 

The  steady  supersonic  two-dimensional  flows  admit  of 
three  families  of  characteristic  lines,  to  wit: 

-► 

--  The  flow  lines  X|i 

—  The  Mach  lines  which  form  with  Xj  the  angles 

*  a  »  4  arc  «tn  ~  ; 


where  a  is  the  velocity  of  sound  in  the  flow  of  velocity  V. 


In  order  to  preserve  a  parallel  with  the  perfect  gases 
by  introducing  a  polytropic  coefficient  of  expansion  IT*  such 
that 


r 


2 

P 


the  thermodynamic  relationship 


dp 

T  d  5  •  d  H  -  -p- 


wlth  H (p,  p )  (enthalpy),  leads  to 


J  H 


r  _ _ 

P  3  H 

p  '  3  Log  p 


If  we  apply  this  result  to  the  equation  of  state  of  a 
solid  of  the  form 

p  ■  p,  *  t  (T  -  T.) 
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In  the  contacting  media  (explosive,  detonation  products, 
metal),  the  flows  cannot  be  treated  separately  since  each  of 
them  reacts  with  that  adjacent  to  it  in  order  to  assure  the 
equality  of  pressures  (p)  and  deflections  (<p)  at  all  points 
situated  on  the  Interface.  Hence  it  is  preferable  to  choose 
these  two  variables  for  the  solution  of  the  problem. 

Then  the  relationships  along  the  characteristics  are 
expressed  by  (Ref.  ?.  Carridre,  course  of  the  CESM*): 


(4 

cl) 

tin 

a  cot 

a  3p 

rp 

dr) 

dT) 

<- 

a) 

sin 

a  cos 

a  _d£ 

4  12. 

rp 

H 

dS 

*  0 

and  the  curvature  of  a 

flow 

line 

is 

given 

by 

k  .  Jj.  .  (llL  .  iE.) 

2Tp  V3f,  dr)  / 


The  solution  Is  carried  out  by  successive  approximations. 

We  have  assumed  that  the  shock  waves  Intervening  in  the 
calcuJations  are  plane,  hence  the  flow  downstream  is  lsentropic. 
Moreover  we  have  assumed  that  the  compression  waves  which  may 
appear  in  the  flow  art  sufficiently  weak  so  as  not  to  lead  to 
shocks,  and  consequently,  to  entropy  variations. 

The  detonation  front  of  the  explosive  Is  assumed  to  be 
plane  and  perpendicular  to  the  interfaces  (explosive  side). 

♦Research  Center  in  Advanced  Mechanics  (Centre  d' Etudes 
Superleures  de  Mecanlque,  rue  P.  Curie,  Paris). 
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Me  neglect  the  reaction  zone  and  assume  that  the  conditions  of 
Chapman- Jouguet  (C.J.)  are  realized  Immediately  behind  this 
front  which  constitutes  the  sonic  line  of  the  flow  of  the  deto¬ 
nation  products.  In  addition  we  assume  that  the  latter  behave 
like  a  perfect  gas  having  a  polytropic  coefficient  of  3. 

While  the  metal  is  compressed,  the  detonation  products 
expand  until  the  equality  of  the  pressures  and  deflections  in 
the  two  media  in  the  vicinity  of  the  interface  is  assured. 

To  this  end  the  metal  is  subjected  to  a  plane  shock  wave, 
making  an  angle  d"  w  th  the  original  direction  of  the  interface, 
while  a  Mayer  expansion,  centered  at  the  point  of  contact  of 
the  letonatlon  front  and  the  interface,  governs  the  flow  of 
the  detonation  products. 


G> 


front  do  ^dtonation 

yrodulfs do 


<D 

axploaif 
sol ido 


ddtonto  do  Mayor  & 
'P  .in  tar  fact 


choc  dans  to  oital\@ 


(1)  -  detonation  front;  (2)  -  solid  explosive;  (3)  - 
detonation  products;  (4)  -  Mayer  expansion;  (5)  -  shock 
in  the  metal;  (6)  -  expansion  of  the  detonation  prod¬ 
ucts;  (7)  -  "heart-shaped"  curve  of  the  metal. 

Let  U  and  u  be,  respectively,  the  velocity  of  the  shock 
and  of  the  material  behind  the  shock  front  in  the  metal  in  a 
reference  system  linked  to  the  laboratory,  v  the  velocity  of 
the  material  in  the  reference  system  linked  to  the  detonation 
front: 


The  various  relationships  linking  the  flow  parameters 

are: 

—  Geometric  relationships  derived  from  the  configuration 
of  the  shock. 

V  •  D  ain  a  — H_  = - 2 -  .  — V_ 

•In  9  cos  (o  -  <(,)  cos  o 

u  -  u  ,  . 

— - —  •  sin  (  -  >  <p) 

—  Equations  of  the  shock: 

p  *  p  Uu  p(U  -  o)  *  i.  v  E  -  F  •  5  p  (/-  -  -) 

0  £  \0o  p  • 

(p0  is  neglected  compared  with  p). 

--  Equation  of  state  of  the  metal. 

f  <f.  p.  T)  --  o 

These  equations  permit,  in  particular,  to  determine  for 
each  value  of  p  the  corresponding  value  of  cf,  when  the  velocity 
of  detonation  D  of  the  explosive  is  given  ("heart-shaped"  curve 
of  the  metal ) . 

The  intersection  of  the  shock  curves  of  the  metal  and 
the  curves  of  the  expansion  of  the  detonation  products  furnishes 
in  plane  (p,  <p)  the  values  which  determine  the  Initial  conditions 
of  flow  in  the  metal  after  shock. 

It  1 s  to  be  noted  that  tensions  (negative  pressures) 
have  appeared  within  the  metal  in  the  course  of  the  study.  They 
were  treated  simply  by  extending  the  calculations  to  valuee  of 
p  greater  than  pQ .  The  calculations,  which  were  carried  out  up 

to  -50  kb,  did  not  reveal  any  anomalies. 

The  numerical  values  employed  during  the  study  of  the 
lifting  of  copper  plate  of  thickness  £  by  ft  plate  of  explosive 
of  thickness  e,  such  that  6/e  =  0.148,  are: 

—  Explosive:  pcj  =  3  x  10?  pi  ezes;  D  =  8,100  m/sec 

—  Metal  (copper):  p0  =  6.92  g/cm^. 
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Coefficients  of  the  P.E.J.  equation  of  the  metal: 

rc=  381  x  105  plezes;  y6  =  IO.831  (nondimensional ) 

g  =  6,809  plezes/OK;  Cv  =  373.3  kj/t 

Initial  conditions  after  shock  in  the  metal: 

=  2  x  107  piezes;  a^  s  4,783  m/sec 

Th  =  403. 5°K;  <oK  =  9.9457  g/cm3 

<r  a  35° 1 5  *  (angle  between  the  shock  and  the  metal/ 
explosive  interface) 

«p  s  2°55'  (angle  of  initial  deflection) 

V^j  =  7,828.33  ra/sec  (velocity  of  the  material  in  the 

reference  system  linked  to  the 
detonation  front). 

The  critical  angle  obtained  in  this  way  is  in  very  good 
agreement  with  the  experiments  (deviation  of  the  order  of  the 
experimental  error):  Plates  2  and  3* 

3.2.  THE  TRANSITION:  SHOCK  -  DETONATION 

3.2.1.  Validity  of  the  Experimental  Method  Employed 

The  experimental  Btudy  of  the  generation  of  detonation 
by  impact  in  a  solid  granular  explosive  has  been  carried  out  by 
the  wedge  method,  described  in  Chapter  2. 

Thus,  the  one-dimensional  character  of  the  phenomenon  is 
realized  within  the  block  of  explosive  while  it  is  probably  not 
realized  in  the  vicinity  of  the  free  surface  of  the  wedge  where 
the  measurement  is  carried  out.  In  effect,  in  this  region  the 
shock  is  not  plane  and  Its  form  changes  during  the  entire 
transitory  period  which  precedes  the  establishment  of  the 
steady  state  of  detonation,  since  the  conditions  downstream  of 
the  shock  front  constantly  charge  during  this  period.  Hence 
the  experimental  measurement  is  fundamentally  distorted  by 
errors. 


A  theoretical  study  of  this  influence  is  impossible  at 
the  present  time  since  the  phenomenon  is  three-dimensional 
(two-dimensional  if  the  medium  is  homogeneous),  unstationary 
with  chemical  kinetics. 


fac*  d'lapact  <i) 


(1)  -  impact  surface;  (2)  -  explosive  wedge; 

(3)  -  perturbed  zone;  (4)  -  one-dimensional 
character  respected. 

Despite  this  serious  defect,  this  method  is  the  only 
one,  to  our  knowledge,  which  permits  a  continuous  recording  of 
the  variation  of  the  reactive  shock  in  the  explosive. 

Although  the  results  are  different  in  absolute  value 
when  the  angle  of  the  wedge  is  varied,  it  is  nevertheless 
feasible  that  for  the  sane  angle  they  should  be  directly 
comparable. 

A  second  difficulty  is  due  to  the  impact  Itself.  Since 
the  experiments  were  carried  out  at  atmospheric  pressure,  some 
air  is  captured  between  tne  projectile  and  the  target.  Hie 
Influence  of  this  gas  13  twofold: 

—  It  acts  like  a  piston  to  assure  that  the  target  is 
gradually  brought  into  motion;  at  the  time  of  the  Impact,  the 
target  no  longer  has  zero  velo>  ity; 

—  This  air  has  a  high  temperature  as  a  result  of  the 
numerous  reflections  of  the  shock  between  the  target  and  the 
projectile.  Consequently,  before  the  pressure  signal  furnished 
by  the  impact,  the  explosive  is  subjected  to  a  high- intensity 
thermal  signal. 

The  influence  of  the  air  which  acts  like  a  piston  con¬ 
tributes  only  a  slight  advance  into  the  determination  of  the 
initial  moment  of  entry  of  the  signal  into  the  target.  This 
error  is  eliminated  once  the  shock  due  to  the  impact  overcomes 
the  compression  waves  Induced  by  the  compressed  air. 
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As  for  the  thermal  pulse,  theoretical  studies  i.  10]  show 
that  It  does  not  have  the  time  to  modify  the  behavior  of  the 
explosive  on  account  of  the  low  thermal  conductivity  of  this 
substance. 

3.2.2.  Practical  Arrangement 

The  experimental  device  comprises  two  parts: 

—  An  argon  chamber  placed  in  the  plane  of  Impact  of  the 
projectile ; 

--  An  explosive  prism  cut  out  of  the  mass,  whose  cross 
section  1 8  an  Isosceles  triangle.  Two  of  its  faces  are  covered 
with  mylar  (10  p  thick),  which  has  been  aluminized  so  as  to 
reflect  the  luminosity  produced  during  the  experiment  by  argon 
flashes . 


C  *  16*45'  ,,4n  j) 


(1)  -  impact  plan;  (2)  -  surfaces  covered  with  aluminized  "my¬ 
lar";  (3)  -  Impact  surface;  (A)  -  argon  flash;  (3)  -  toward 
camera. 
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This  device,  which  permits  the  simultaneous  examination 
of  the  shock  on  two  faces  of  the  explosive  prism,  makes  it 
possible  to  take  into  account,  during  the  treatment  of  the 
results,  the  slight  obliqueness  of  the  projectile  with  respect 
to  the  impact  plane. 

The  recording  is  made  by  means  of  a  slit  camera. 

The  whole  assembly  of  the  projection  device  plus  target 
is  shown  on  Plate  4. 

It  is  found  experimentally  that  when  the  shock  attains 
a  point  of  the  free  surface  of  the  explosive,  the  deflection  of 
the  mylar  stops  abruptly  the  reflection  of  the  luminosity 
originating  from  the  argon  flash.  This  is  reflected,  for  the 
whole  of  the  target  explored  by  the  slit  of  the  camera,  by  a 
curve  which  is  detached  in  black  on  a  white  background  given 
before  the  passage  of  the  reactive  shock  by  the  argon  flash, 
and  after  its  passage  by  the  detonation  products. 


(1)  -  progress  diagram  of  the  reactive  shock; 

(2)  luminosity  due  to  the  argon  flash;  (3)  - 
zone  of  intense  luminosity  due  to  the  detona¬ 
tion  products. 

Such  an  assembly  makes  it  necessary  to  place  in  the 
image  plane  of  the  input  objective  of  the  slit  camera  not  one 
slit  as  In  the  usual  case  but  two,  the  first  for  the  argon 
chamber  and  the  second  for  the  explosive  prism. 

Hence  on  the  photographic  plate,  one  obtains: 


—  The  visualization  of  the  planeness  of  the  projectile 
on  impact  as  well  as  the  angle  at  which  it  strikes  the  target; 

—  The  progress  diagram  of  the  gradually  reactive  shock 
front,  which  is  propagated  on  the  two  sides  of  the  explosive 
prism; 


—  The  magnification  of  the  optical  assembly,  obtained 
by  preliminary  recording  of  two  luminous  points  placed  in  the 
impact  plane  of  the  projectile. 

Finally,  various  marks  placed  on  the  argon  chamber  make 
it  possible  to  determine  the  extremities  of  the  explosive  prism 
with  precision. 

3.2.3*  Control  of  the  Conditions  of  Impact  on  an  Inert  Target 

To  verify  whether  the  experimental  device  is  satisfac¬ 
tory,  we  have  carried  out  experiments  in  which  the  explosive 
constituting  the  target  was  replaced  by  an  inert  substance  — 
copper  or  aluminum  (AU4G). 

When  the  projectile  and  the  target  are  of  the  same  sub¬ 
stance,  the  determination  of  the  conditions  of  shock  is  im¬ 
mediate.  In  effect,  the  measurement  of  the  impact  velocity 
(Vp)  furnishes  the  velocity  of  the  material  u,  u  =  Vp/2,  while 

the  pressure  is  obtained  by  means  of  p  =  p0  U  u,  where  U  Is  the 

shock  velocity  determined  experimentally. 

These  calculations  are  valid  if  the  projectile  and  the 
target  are  under  the  same  thermodynamic  conditions  at  the 
moment  of  impact.  The  projection  procedure  used  probably  does 
not  ensure  this  possibility,  and  the  projectile,  even  if  it  is 
at  zero  pressure  at  the  moment  of  Impact,  Is  probably  not  at 
the  same  temperature  as  the  target. 

A  theoretical  study  was  undertaken  in  order  to  take 
this  difference  into  account,  and  the  projectile  was  assumed 
to  be  at  a  temperature  of  500°K  instead  of  300°K: 

—  The  determination  of  the  shock  polar  of  the  copper 
from  the  initial  conditions  (p0  =0;  TQ  =  500OK;  =  8.842 

g/cm3  instead  of  8.92  g/cm3  for  300°K)  while  preserving  the 
same  coefficients  for  the  equation  of  state  shows  that  the 
pressure  Is  1  to  2%  lower  than  that  In  the  case  of  300°K  for 
the  same  value  of  material  velocity; 
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--  Calculation  or  this  polar  from  the  coefficient  of 
expansion  of  copper  furnishes  a  much  lower  deviation,  of  the 
order  of  0. 1 — 0.2  «. 

Since  the  experimental  dispersion  Is  greater  than  these 
values,  It  is  not  possible  to  Judge  wht*  ler  it  Is  necessary  to 
take  into  account  the  temperature  rise  of  the  projectile  during 
the  time  that  it  is  set  to  motion  in  order  to  calculate  the 
conditions  of  the  shock  (see  diagram  6,  a  detailed  example  of 
the  treatment  of  the  experiment.!  data  with  the  three  shock 
polars,  to  calculate  the  shock  conditions). 

When  the  projectile  and  the  target  are  made  of  two  dif¬ 
ferent  materials,  the  proalem  is  more  complicated.  It  is  in 
effect  necessary  to  know  the  shock  polar  of  the  projectile  to 
determine  the  conditions  of  impact  In  the  target. 

The  small  deviation  observed  above  between  the  different 
shock  polars  calculated  for  copper  makes  it  possible  for  us  to 
carry  out  the  data  treatment  with  only  one  of  them.  We  chose 
that  which  verifies  the  conditions:  pn  =  0;  Tj  =  300°K. 

Moreover  the  projectile  does  not  simultaneously  strike 
the  entire  target  surface,  3lnce  an~le  nt of  the  experimental 
setup  is  not  strictly  equal  to  y ,  the  angle  of  lifting  of  the 
plate.  We  have  taken  this  slight  obliqueness  on  Impact  into 
account  and  determined  V  by  the  equation 

-2  1)  sin  *  ros  (  j  ^  ) 

Instead  of  Vp  -  D  sin  which  is  valid  only  when  ql  =  <p. 

The  velocity  of  the  shoex  (U)  In  the  target  was  deter¬ 
mined  In  two  ways. 

--  3y  means  of  the  time  which  separates  the  Impact  of 
the  projectile  on  the  two  ends  of  the  target  (which  may  be  used 
when  rc  Is  sufficiently  different  from  f); 

X  r 

r  tg  Y  *  pp-j;  tg  a 

V  =  — - 

1 81“  (J)  V(1  -  r):  +  r!  tgJ  ' 


(for  notations,  see  diagram  below). 


—  By  means  of  the  registering  of  the  diagram  of  prog¬ 
ress  of  the  shock  in  the  target: 

,  V  «1«  (P  *  Y) 
k  tg  9  cos  p 


where  V  is  the  sweeping  velocity  of  the  camera,  and  k  the  mag¬ 
nification  of  the  installation. 


K««ure  do  It  Vitesse 
d'iapact  et  controls 
de  plane  its  du 
jectile  Q 


pro- 


En reaistrement 
du  choc  dans 
la  cibla 


(1)  -  photographic  recording;  (2)  -  measurement  of 
the  velocity  of  impact,  and  control  of  the  planeness 
of  the  projectile;  (3)  -  recording  of  the  shock  in 
the  target;  (4)  -  section  of  target. 

Such  a  type  of  recording  is  shown  in  Plate  5»  Pig-  1. 

The  results  furnished  by  the  two  methods  are  in  good 
agreement,  and  the  values  of  p,  u  obtained  are  shown  in  diagram 
7  and  Table  8. 

It  is  noted  that  the  conditions  of  impact  determined  in 
this  way  are  situated  in  a  rather  satisfactory  manner  on  the 
shock  polars  derived  by  other  experimental  methods  (cf.  L7l, 
for  example). 
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3. 2. 4.  Experimental  Results  Obtained  with  Explosive  D 

We  have  seen  that,  In  order  to  Initiate  the  detonation, 
the  Intensity  and  duration  of  application  of  the  pressure  sig¬ 
nal  induced  In  the  target  explosive  should  he  greater  than  the 
minimum  values  pm  and  Z . 

In  the  experimental  study  which  follows,  only  the  first 
point  wa9  examined,  even  though  the  device  selected  permits 
varying  the  two  parameters.  However,  our  current  knowledge  of 
the  behavior  of  metals  at  relatively  low  r"-c3sures  (below  100 
kb)  do  not  permit  relating  the  duration  of  the  signal  Induced 
In  the  target  to  the  thickness  of  the  projectile.  In  effect, 
the  assumption  which  we  have  made  --  of  a  behavior  governed  by 
the  laws  of  hydrodynamics  (perfect  fluid)  —  Is  no  longer 
acceptable  on  account  of  the  non-negliglble  Influence  of  the 
elastoplastlo  character  of  the  metal  constituting  the  projec¬ 
tile.  The  present  deviation  between  hydrodynamic  theory  and 
experiment  Is  very  large.  In  certain  cases,  the  duration  of 
application  of  the  pressure  signal  Is  but  one  half  of  the 
theoretical  time  (cf.  Chapter  5,  [6]). 

In  order  to  vary  the  impact  velocity  (Vp),  we  changed 
the  thickness  (<£)  of  the  projectile  while  maintaining  the  char¬ 
acteristics  of  the  projecting  explosive  (material  and  thick¬ 
ness)  constant. 

This  procedure  has  permitted  us  to  obtain  the  following 
Impact  velocities  with  a  copper  projectile,  while  the  theoret¬ 
ical  determination  of  the  shock  polar  of  the  Inert  explosive 
from  the  experimental  results  (cf.  Chapters  1  and  5)  has  fur¬ 
nished  us  the  Intensity  of  the  shock  Induced  in  the  target 
(diagram  9). 

The  values  obtained  are  listed  In  the  table  below: 


Epaisseur  du  projectile  i 

(mm)  j) 

2 

J _ 

3 

4 

5 

Vitesse  d'impact  Vjim/s) 
(a  1  30  m/s) 

1  870 

1500 

1200 

1000 

Pression  induite  (kb) 

(S' 

131 

96 

70,5 

55,5 

(1)  -  Thickness  of  projectile;  (2)  -  Impact  velocity 
(to  within  +  30  ra/sec);  (3)  -  Induced  pressure. 
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The  precision  with  respect  to  the  preBeure  la  certainly 
quite  low  due  to  the  experimental  aiificulties  encountered 
during  the  determination  of  the  dynamic  adiabatic  of  the  Inert 
explosive. 

Experiments  carried  out  at  800  ra/sec  have  not  permitted 
us  to  observe  tho  detonation  of  explosive  D  on  impact,  since 
no  luminosity  has  been  recorded  by  the  photographic  plate. 

The  results  obtained  for  impact  velocities  greater  than 
this  value  show  that  the  latter  is  a  fundamental  parameter  of 
the  phenomenon  under  investigation. 

While  for  Vp  =  1,000  m/sec  the  steady  detonation  state 

is  attained  In  approximately  3  fisec,  this  time  is  reduced  to 
1  psec  when  Vp  =  1,200  m/sec.  This  value  then  varies  very 
little  when  the  velocity  of  impact  Increases,  as  can  be  seen 
on  the  photographic  plates,  since  the  zone  in  question  Is  then 
very  small  and  the  smallest  deviations  of  the  planeness  of  the 
projectile  notably  perturb  the  recording. 

It  should  be  remarked  that  In  this  type  of  experiment 
the  moment  of  Impact  of  the  projectile  on  the  target  is  not 
defined  with  precision.  It  cannot  be  given  by  the  argon 
chamber  which  has  a  certain  "response  time,"  and  the  start  of 
the  diagram  of  the  progress  of  the  shock  in  the  t  irget  is  dif¬ 
ficult  to  read  off,  as  we  shall  see  be’ow. 

With  regard  to  all  these  recordings  we  can  make  the 
following  observations  which  are  more  obvious  as  the  impact 
velocity  is  lower: 

--  Before  Impact  the  extreme  edges  of  the  target  no 
longer  reflect  the  light  furnished  by  the  argon  flash.  This 
may  be  due  to  the  presence  of  the  air  cushion  which  precedes 
the  projectile.  This  air  cushion  may  either  perturb  the 
target  explosive  by  the  generation  of  compression  waves  in  the 
latter's  interior,  or  project  beyond  the  side  of  the  target 
and  deviate  the  light  beam  originating  from  the  argon  flash. 
Nevertheless,  the  thickness  of  the  explosive  subject  to  this 
effect  is  small,  of  the  order  of  a  millimeter,  and  the  length 
of  time  during  which  this  perturbation  makes  itself  felt  is 
of  the  order  of  2  |usec.  Under  these  conditions  it  is  impos¬ 
sible  to  define  precisely  the  moment  of  Impact  of  the  projec¬ 
tile  on  the  target; 

--  Immediately  upon  Impact  the  luminosity  due  to  the 
argon  flash  is  abruptly  stopped,  and  the  phenomenon  Is  analogous 
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to  the  case  of  an  Inert  target,  despite  the  fact  that  the 
progress  diagram  recorded  shows  a  slight  acceleration  of  the 
shock; 


—  Finally,  the  explosive  becomes  luminous  as  soon  as 
the  shock  wave  arrives;  a  distinct  slope  discontinuity  appears 
on  the  recording  and  the  velocity  of  displacement  of  the 
recorded  phenomenon  gradually  approaches  the  value  of  the 
steady  detonation  state. 

For  reasons  mentioned  above  the  moment  of  Impact  is  not 
known  accurately.  We  have  determined  it  by  extrapolation  of 
the  curve  representative  of  the  shock  on  the  photographic  plate, 
in  ordor  to  compare  the  various  recordings  of  the  same  series 
of  experiments.  Then  it  may  be  thought  that  for  two  Identical 
assemblies,  this  moment  is  defined  to  within  a  constant,  while 
the  time  deviations  shown  or  the  diagrams  are  determined  with 
the  following  precision:  the  deviation  between  two  points  of 
the  recording  is  measured  to  within  ♦  5/ 100  mo,  which  cor¬ 
responds  to  a  determination,  in  time,  to  within  ♦  10  nanoseconds 
(sweeping  velocity  6.4  mn/^aec)  and  in  space  to  within  ♦  2/10 
mm  (magnification  of  the  setup  l/*+).  “ 

Some  experiments  have  been  carried  out  with  larger 
explosive  wedges  in  order  to  follow  the  reactive  shock  in  the 
explosive  for  a  longer  period  of  time  (see  Plate  5,  Fig.  3). 

The  results  obtained  in  this  wav  are  shf wn  in  diagrams 
10,  11,  12  and  13. 
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Plate  4 


Assembly  of  the  shock-generating  device  and  target, 
linear  generator;  (2)  -  argon  chamber;  (3)  -  target. 


a  Vp  1.873  in /s  b  :  V,  »  1.208  m/s 


Fig.  1.  Impact  on  an  Inert  target, 
a)  made  of  copper;  b)  made  of  aluminum  (AU  AG). 


double  wedge  single  wedge 

V,  =  1.020  m/s  V,  *  1.000  m/s 


Fig.  2.  Impact  on  a  target  made  of  explosive  D. 

Note:  On  each  photo  the  distance  between  the  two  vertical  linee 
is  200  mm. 
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Diagram  6 
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different  froac  U(0);  (6)  -  Al. 


Diagram 


Diagram  12 
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Chapter  4 

CRITICAL  EXAMINATION  0?  THE  VARIOUS  THEORETICAL  STUDIES  OF 
THE  INITIATION  OF  DETONATION  BY  SHOCK  WAVE 


The  theoretical  study  of  the  Initiation  of  detonation 
by  Impact  or  by  frontal  shocX  transmitted  by  means  of  a  barrier 
Is  essentially  the  same.  We  have  seen  that  the  experimental 
setup  differs  only  In  regard  to  the  shape  of  the  pressure  sig¬ 
nal  Induced  In  the  explosive.  The  signal,  furnished  by  the 
Impact  of  the  projectile  In  the  first  case,  Is  of  rectangular 
shape  at  a  given  Instant  In  plane  x,  p,  while  the  signal  trans¬ 
mitted  by  the  Inert  barrier  In  the  second  case  Is  of  triangular 
shape. 


The  present  theoretical  schemes  do  not  permit  giving  an 
account  of  the  behavior  of  the  heterogeneous  3olid  explosive, 
whereas  In  the  case  of  homogeneous  explosives  they  furnish 
results  which  are  In  rather  good  agreement  with  experiments. 

In  a  system  using  the  Lagrange  variables  (X,  t),  the 
number  of  functions  which  we  must  determine  Is  seven  In  the 
most  simple  case:  u,  x,  p  ,  E,  p,  T  and  m  (for  notations,  see 
table  at  the  end  of  the  chapter). 

The  four  hydrodynamic  equations  must  be  completed  by 
three  other  equations  In  order  to  be  able  to  solve  completely 
the  system  of  seven  functions  with  two  variables. 

Two  of  these  three  equations  Involve  the  behavior  of  the 
medium;  of  the  two,  one  Is  an  equation  of  state  f(p,  v,  T)  r  0, 
while  the  other  permits  the  determination  of  Its  Internal 
energy  for  all  values  of  (p,  v,  T),  or  E(p,  v,  T)  =  0. 

The  last  relationship  Is  the  law  of  liberation  of  chem¬ 
ical  energy,  which  furnishes  the  reaction  rate  of  the  explosive. 
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4.1.  HYDRODYNAMIC  EQUATIONS 

In  a  system  of  Lagrange  coordinates,  the  three  flow 
equations  are  written  as  follows,  if  the  viscosity  and  thermal 
conduction  are  neglected: 


ax 

TT 


P. 


3u  .  3p  „ 

p,  TT  *ax  *  0 


it — TT  4  P1T  0 


(conservation  of 
(conservation  of 
(conservation  of 


mass ) 
momentum ) 
energy ) 


To  these  three  equations  we  have  to  add  the  kinetic 
equation: 


U 


which  links  the  particle  velocity  to  Euler's  variable 
LX ( X,  t ) ] . 

The  phenomenon  which  we  want  to  study  has  a  duration  of 
the  order  of  a  few  microseconds,  hence  neglecting  the  viscosity 
and  the  heat  conduction  is  a  Justified  assumption,  since  the 
effects  for  which  they  are  responsible  are  much  slower  than 
those  produced  by  the  dynamic  parameters  (a  few  milliseconds 
Instead  of  a  few  microseconds). 

Nevertheless,  Enig  takes  these  two  factors  into  account 
in  his  theoretical  study,  .'hlle  expressing  some  doubt  as  to 
their  physical  meaning.  Their  Interest  lies  rather  in  the 
Interpretation  which  they  give  to  Rlchtmeyer's  concept  of 
pseudoviscosity ,  Introduced  in  the  system  of  hydrodynamic 
equations  for  the  treatment  of  shocks  [4]. 

4.2.  EQUATION  OF  STATE 

The  simplest  system,  containing  seven  equations, 
requires  the  use  of  a  single  equation  of  state  to  represent 
the  solid  explosive,  the  detonation  products  and  the  gradual 
passage  from  one  to  the  other. 

Let  us  see  what  knowledge  we  have  of  these  various 

states . 
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The  dynamic  methods  now  make  It  possible  to  know,  with 
good  precision,  the  dynamic  adiabatic  of  an  Inert  solid,  and 
to  derive  from  It  the  solid's  equation  of  state,  at  least  In  a 
limited  area  of  Clapeyron's  plane.  However  we  are  reduced  to 
guesses  when  this  substance  Is  capable  of  "reacting."  In 
effect,  as  we  have  seen  above,  the  measures  carried  out  In 
this  case  are  then  distorted  by  the  start  of  the  chemical  reac¬ 
tion.  Only  those  results  which  have  been  obtained  at  low  pres¬ 
sure  may  be  considered  valid. 

For  the  detonation  products,  the  equation  of  state  Is 
known  sufficiently  well  only  at  the  Chapman- Jouguet  (C.J.) 
point.  It  should  be  noted,  moreover,  that  the  study  of  the 
ieentroplc  expansion,  which  takes  into  account,  for  each 
pressure,  the  new  conditions  of  thermodynamic  equilibrium  of 
the  constituent  gases,  Is  in  rather  good  agreement  with  the 
experiments . 

Finally,  there  are  no  uata  which  would  make  it  possible 
to  form  an  Idea  of  the  behavior  of  the  solid-gas  mixture  during 
the  reaction.  Nevertheless  we  assume  that  this  reaction  zone 
may  be  divided  Into  small  regions  In  which  thermodynamic  equil¬ 
ibrium  Is  established,  making  the  unequivocal  definition  of  the 
parameters  of  state  possible  In  each  medium  (pressure,  specific 
volume,  temperature). 

Hubbard  and  Johnson  [7]  --  the  first  researchers  to  have 
proposed  a  mathematical  model  on  this  subject  --  employ  a 
single  equation  of  state  valid  for  the  solid  and  for  the  deto¬ 
nation  products: 

p(v  -  b)  •  n  K  T 


This  equation  Is  that  of  a  gas  whose  internal  energy  Is 
a  function  only  of  the  temperature:  E  =  CyT.  This  eystem 
obviously  does  not  take  Into  account  the  state  of  advance  of 
the  reaction  and  can  only  furnish  a  highly  approximate  agree¬ 
ment  with  experiments. 

In  particular,  although  this  equation  of  state  repre¬ 
sents  the  detonation  products  rather  well,  It  furnishes  for' 
the  Initial  conditions  of  the  solid  explosive  a  particularly 
low  temperature,  not  related  to  reality.  In  effect,  the 
equation  p(v  -  b)  =  nRT  must  be  satisfied  both  by  the  solid 
explosive  before  the  reaction  pn ,  vn ,  Tn ,  and  by  the  detonation 

^  ^  ^ 

products  (C.J.  conditions:  p,  v,  T). 
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Hence  we  rauet  have,  with  b  and  nR  being  constant, 


R  (v.  -  *»)  -  it  .1  T.  A  -  b)  •  n  R  f 


or 


■  b) 

TT 


How,  the  denominator  of  the  right  hand  side  of  this 
expression  Is  very  much  greater  than  the  numerator,  since  the 

vn  -  b  A 

numerical  value  of  -  Is  only  a  few  units  while  p  Is  of 

v  -  b 

the  order  of  200,000  pn  Li 2]. 

A  numerical  application,  using  the  values  associated 
with  the  explosive  employed  In  the  calculations  below,  gives 
the  following  results: 

T,  •  0.057*  K 

b  •  0.430  v.  •  0,588  0  -  0.480  (cm>/«) 

P  •  216  kb  A  t  •  4.124*  K 


The  use  of  such  an  equation  entails  serious  difficulties 
for  the  Initiation  of  the  reaction  if  we  assume  that  this 
Initiation  Is  a  function  of  the  temperature  through  the  Inter¬ 
mediary  of  Arrhenius'  exponential  law. 

Hence  It  Is  desirable  to  look  for  a  way  of  improving 
the  representation  of  the  transition:  solid  — >  gas.  The  use  of 
a  single  equation  --  even  If  It  takes  Into  account  the  degree 
of  advance  of  the  reaction  --  cannot  lead  to  satisfactory 
results.  Hence  It  Is  preferable  to  look  for  a  model  which 
employs  one  equation  of  state  for  each  of  the  two  media,  the 
solid  and  the  gas.  This  method,  however,  considerably  compli¬ 
cates  the  mathematical  model. 

Assuming  that  there  Is  a  single  reaction  having  the 

form 


Solid  explosive  — »  Detonation  products 

with  an  equation  of  state  for  each  medium,  we  Introduce  the 
variables  pe ,  vg,  Tg ,  Eg  and  p?,  vf,  Tg,  Ep,  with  subscript  e 

denoting  the  solid  medium  and  g  the  gaseous  medium. 
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To  simplify  the  mathematical  model,  various  authors  have 
assumed  in  their  study  of  homogeneous  explosives  that  during 
the  entire  reaction  pe  =  pg  =  p  and  T8  =  Tg  =  T  [4,  10,  1 1 3 « 

Let  us  s'e  what  these  simplifications  represent,  and  examine 
whether  it  is  possible  to  accept  them  for  the  study  of  hetero¬ 
geneous  granular  explosives. 

At  equal  pressures  (of  the  order  of  the  C.J.  pressure), 
the  temperature  of  a  gas  is  much  higher  than  that  of  a  solid 
(for  example  4,000°K  instead  of  1,800°K)  —  see  the  numerical 
application  to  the  explosive  used  in  the  calculations  below). 

The  equality  of  temperatures  can  only  be  conceived  if  the 
thermal  conduction  permits  the  establishment  of  equilibrium 
between  the  two  media  within  a  short  time  (meaning  short  with 
respect  to  a  microsecond).  Now,  we  know  that  this  is  not  so 
at  all  l5],  since  the  explosive  has  hardly  the  time  to  heat  up 
during  the  few  microseconds  during  which  the  initiation  of 
detonation  takes  place.  Hence  the  assumption  Ts  =  Tg  cannot  be 

accepted  for  the  study  of  the  behavior  of  granular  explosives, 
while  it  is  fully  Justified  in  the  case  of  homogeneous 
explosives . 

The  second  simplification,  equality  of  pressures,  while 
still  debatable,  is  nevertheless  closer  to  reality.  In  effect, 
this  equality  is  linked  to  the  wave  velocities  in  the  two  media: 
solid  and  gaseouB.  For  neither  of  these  media  is  the  velocity 
of  sound  measurable,  and  it  is  necessary  to  make  some  highly 
debatable  assumptions  to  determine  it  theoretically. 

For  the  gas,  whose  pressure  constantly  changes  during 
the  reaction,  a  value  of  a  few  millimeters  per  microsecond 
seems  to  be  of  a  good  order  of  magnitude,  if  we  refer  to  the 
sole  known  value,  that  obtained  theoretically  in  the  Chapman- 
Jouget  plane  (D  -  u  =  a).  In  the  numerical  application  of  the 
next  chapter  we  obtain  approximately  6.7  mm/jisec. 

In  the  case  of  the  solid,  let  us  remember  that  the 
grains  are  of  different  chemical  composition.  Hence  the  deter¬ 
mination  of  a  mean  value  is  a  delicate  undertaking.  3y  anal¬ 
ogy  with  homogeneous  solids  it  may  be  estimated  that  it,  too, 
is  of  the  order  of  a  few  millimeters  per  microsecond. 

The  velocity  of  sound  in  each  of  the  media  may  therefore 
be  estimated  to  be  of  the  same  order  of  magnitude.  Under 
these  conditions,  if  we  assume  that  a  Dressure  equilibrium  may 
be  established  for  each  of  them  in  a  negligible  time  Jn  each 
of  the  regions  which  constitute  the  space  where  the  reaction 
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takes  place,  we  are  Justified  in  likewise  assuming  that  the 
two  media  —  solid  and  gas  —  are  in  pressure  equilibrium. 

Hence  the  equality  p8  =  pg  is  acceptable  in  the  case  of  granu¬ 
lar  explosives. 

The  introduction  of  the  new  functions  p,  v,  T,  E,  of  the 
solid  (s)  and  of  the  detonation  products  (g)  requires,  during 
the  moments  when  these  two  media  are  in  contact,  the  determin¬ 
ation  of  the  mean  values  of  v  and  E  which  enter  into  the  hydro- 
dynamic  equations. 

These  values  are  defined  by 


v  ■  (1  -  m)  v,  4  m  v, 
E  •  (1  -  m)  E,  4  m  E, 


Enlg  [4]  has  selected,  as  equations  of  state,  Talt's 
equation  generalized  for  the  solid  (p  ♦  B)  v  -  (pQ  ♦  b)  v0  = 

a  (y  -  1)  (E  -  )  and  an  equation  of  polytropic  gas  (y  a 

=  constant  for  the  detonation  products.  Mader  has  improved  the 
preceding  scheme  by  using  equations  of  state  that  are  better 
adapted:  Crunelsen's  equation  lor  the  solid,  the  equation  of 
Flckett  and  Wood  for  the  detonation  products.  These  equations 
have  allowed  him  to  obtain,  in  the  case  of  homogeneous  explo¬ 
sives,  an  excellent  agreement  between  theory  and  experiment 
[10]. 

4.3.  CHEMICAL  KINETICS 

The  last  equation  determines  the  rate  of  liberation  of 
the  chemical  energy  of  the  explosive.  This  is  in  a  way  its 
rate  of  reaction.  The  solid  explosive  and  the  detonation 
products  are  treated  as  homogeneous  bodies,  and  it  is  assumed 
that  the  reaction  goes  to  completion. 

Solid  explosive  — >  detonation  products 

The  law  of  chemical  kinetics  most  frequently  used  is 
that  of  Arrhenius  of  the  type: 

~  ■  v  (i  -  m)  • -tr¬ 
ibe  use  of  this  law  is  based  essentially  on  the  fact 
that  we  are  in  the  presence  of  a  chemical  reaction.  The  ac¬ 
tivation  energy  used  is  obtained  by  extrapolation  to  the 
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detonation  temperatures  from  measurements  of  the  rate  of  de¬ 
composition  of  explosives,  carried  nut  at  ambient  temperature 
(at  the  most  a  few  hundred  degrees  centigrade).  The  measure¬ 
ments  carried  out  by  various  laboratories  on  a  large  number  of 
explosives  lead  to  an  activation  energy  which  Is  not  too 
variable  --  30  to  50  kcal/mole  --  but  they  lead  to  frequency 

factor  (v)  situated  between  lO1^  and  10^0  see-'  [5]. 

This  extrapolation  to  the  detonation  temperatures  of  the 
results  obtained  at  a  few  hundred  degrees  Is  highly  debatable, 
since  the  chemical  phenomenon  is  not  governed  strictly  by  a 
single  equilibrium  equation  but  by  a  number  of  equations  such 

as 


co,  *■  c  2  co 
CO,  +  H,  ; - »co  -*  h,o 

2  H,0  , - »  2  H,  ♦  O, 


Of  all  these  reactions  the  slowest  one  Is  that  which, 
at  a  given  Instant,  Imposes  Its  velocity  as  a  function  of  the 
thermodynamic  conditions  of  the  medium.  However  these  con¬ 
ditions  constantly  change  as  a  result  of  the  evolution  of  the 
"reaction,"  and  there  Is  no  assurance  that  the  slowest  of  them 
remains  the  slowest  fer>  a  given  pressure  and  temperature 
regardless  of  the  magnitude  of  these  two  parameters. 

On  the  other  hand,  tne  frequency  factor  which  may  be 
obtained  from  these  experiments  varies  to  a  considerable  extent 

(♦  10^).  This  then  makes  It  possible  to  adjust  the  theoretical 
and  experimental  results  in  a  relatively  easy  manner.  Then 
the  conclusion  drawn  from  this  that  the  order  of  magnitude  of 
the  "Incubation  time"  of  homogeneous  explosives  found  by  cal¬ 
culation  Is  compatible  with  that  obtained  experimentally  is 
not  a  surprising  one. 

Because  of  its  exponential  term,  this  law  Is  particularly 
sensitive  to  the  temperature  and  if  we  do  not  want  to  have  an 
abrupt  Increase  of  the  reaction  rate  for  a  small  temperature 
Increase,  we  must  use  In  the  numerical  calculations  very  small 
A*  and  At  steps.  Choosing  a  Ax  of  the  order  of  a  micron  for 
a  homogeneous  substance  is  acceptable,  but  this  value  seems  to 
be  low  for  a  heterogeneous  explosive  whose  averace  grain  size 
Isa  hundred  times  greater. 
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(1)  -  Temperature  in  the  detonation 
wave;  (2)  -  extrapolation  of  the  reac¬ 
tion  time  to  the  detonation  temperatures; 

(3)  -  [5,  P  152];  (4)  -  experimental 
region. 

The  exponential  character  of  this  law  is  reflected  also 
by  the  fact  that  a  temperature  difference  of  a  few  tens  of 
degrees  is  sufficient  to  bring  about  or  prevent  the  initiation 
of  the  reaction.  This  temperature  variation  may  be  obtained  in 
two  ways,  by  modifying 

—  The  intensity  of  the  initiating  shock-wave, 

—  The  initial  temperature  of  the  explosive. 

While  this  law  explains  the  behavior  of  homogeneous 
explosives  in  a  satisfactory  manner  —  in  this  case  it  is  found 
that  above  a  certain  pressure-  (or  temperature-)  threshold  the 
reaction  is  violent  --  it  does  not  make  it  possible  to  account 
for  the  much  more  gradual  variation  of  the  initiation  of  deto¬ 
nation  in  a  granular  explosive  when  the  temperature  generated 
by  the  shock  wave  is  modified. 

Finally,  the  exponential  term  of  Arrhenius'  equation 
Intervenes  for  the  determination  of  the  initial  point  of  the 
complete  reaction.  In  effect,  after  the  passage  of  the  shock, 
the  amount  of  chemical  energy  liberated  is  very  small.  Only, 
through  the  cumulative  effect  which  it  gives  rise  to,  time 
ensures  the  continual  increase  of  the  liberated  energy.  It 
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can  be  seen  then  that  It  Is  the  Interface,  the  region  which  is 
subjected  to  this  heating  for  the  longest  period  of  time,  that 
attains  the  temperature  above  which  the  abrupt  character  of 
Arrhenius'  law  makes  itself  felt.  This  particularity  is  so 
sudden  that  it  does  not  give  the  hydrodynamic  phenomena  time  to 
attenuate,  in  the  explosive,  this  abrupt  pressure  Increase 
which  would  be  slowed  down  by  a  beam  of  expansion  waves  while 
compression-  or  shock  waves  would  be  propagated  in  the  shock 
generator. 

Nevertheless  it  is  found  experimentally  in  the  case  of 
the  d<  vice  with  barrier  that,  since  the  pressure  signal  is  not 
constant  as  a  function  of  the  time,  the  liberation  of  chemical 
energy  ought  to  slow  down  at  the  interface.  The  measurements 
which  have  been  carried  out  L 3 ]  show,  however,  that  in  the 
case  of  the  homogeneous  explosives  the  detonation  Indeed 
begins  in  this  area.  Consequently  the  decrease  of  pressure, 
and  hence  that  of  temperature,  is  not  sufficient  to  compensate 
for  the  effect  of  time  on  the  liberation  of  the  chemical 
energy  of  the  explosive. 

All  these  remarks  lead  us  to  prefer  a  more  gradual  law 
of  energy  liberation,  not  of  an  exponential  character,  in  our 
attempt  to  explain  the  behavior  of  granular  explosives. 

For  the  sake  of  convenience  we  may,  as  suggested  by 
Klstiakowsky  i9]  choose  a  law  of  energy  liberation  which  is  a 
function  of  the  pressure,  a  parameter  which  intervenes  explicitly 
in  the  hydrodynamic  equations,  instead  of  linking  it  to  the 
temperature. 

This  choice  is  quite  arbitrary,  since  these  two  variables 
are  not  Independent.  Their  direction  of  variation  is  the  same, 
and  during  the  reaction,  the  knowledge  of  one  of  them  determines 
the  other. 

This  idea  haB  been  recently  taken  up  by  Adams  [ 1 ] ,  and 
numerical  calculations  have  been  carried  out  by  V.arner  [l2]  by 
adding  a  pressure  term  to  Arrhenius'  law.  Then  the  law  of 
liberation  of  chemical  energy  has  the  following  form: 


It  makes  it  possible  to  obtain  a  pressure  rise  as  soon  as  the 
shock  passes,  in  contrast  with  the  results  obtained  with 
Arrhenius'  law  alon®,  but  it  does  not  assure  the  progressive 
variation  of  the  "Induction”  time  as  a  function  of  the  intensity 


£0 


of  the  pressure  signal.  The  reason  for  this  is,  probably,  that 
the  exponential  terra  becomes  predominant  during  the  few  moments 
which  precede  the  detonation,  whereas  it  is  the  pressure  term 
which  Imposes  the  rate  of  reaction  during  the  first  part  of  the 
phenomenon.  The  boundary  of  these  two  influences  should  be 
quite  sharp  and  independent  of  the  initial  conditions.  These 
calculations  have  been  slightly  improved  [Warner,  Discussion 
at  the  9th  Sympotlum  on  Combustion,  p  527 J  by  the  use  of  a  law 
of  energy  liberation  which  is  a  function  of  the  state  of  advance 
of  the  reaction.  Ihis  makes  it  possible  to  find  the  ’’overshoot" 
obtained  experimentally  by  Jacobs  [8]  with  a  particular  explosive, 
compressed  TNT,  without,  however,  granting  to  the  pressure  peak 
obtained  in  this  way  the  certainty  of  being  able  to  vary  as  a 
function  of  the  Intensity  of  the  shock. 

Finally,  the  influence  of  the  charge  density,  insofar 
as  it  decreases  the  "Induction"  delay  when  the  porosity  of  the 
explosive  is  Increased,  does  not  lead  to  the  appearance  of  a 
true  variation  of  these  delays  as  suggested  by  Arrhenius’  law. 


J-T  (B 


Velocity  of  the  shock  as  a 
function  of  the  dietance, 
for  various  densities. 
(1)  -  compressed  TNT  lb  = 
sec];  (2)  -  density. 


Pressure  peak.  Holebuming. 

(3)  -  pressure  (atmospheres); 

(4)  -  time. 


To  conclude  we  shall  mention  that  the  agreement  between 
theory  and  experiment  which  i?  attained  in  certain  specific 
cases  is  no  longer  present  when  the  initial  conditions  of  the 
shock  are  modified,  hence  a  more  gradual  law  of  energy  liber¬ 
ation,  l.e.,  one  without  an  exponential  term  containing  e~I£ 

T 
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or  e~-  (which  la  essentially  the  same  thing)  l2,  6,  12]  seems 
P 

preferable.  In  the  mathematical  model  which  we  are  proposing 
the  energy  is  liberated  linearly  as  a  function  of  the  pressure. 

NOTATION 

x  Euler's  variable  (abscissa  of  grid  X  at  time  t) 

X  Lagrange's  variable  (initial  abscissa  of  a  grid) 
t  time 

p  pressure 

v0  Initial  specific  volume 
v  specific  volume  at  time  t 

T  temperature 

u  material  velocity 

E  Internal  energy 

E#  activation  energy 

Qr  heat  of  reaction 

ra  mass  of  gas  formed  by  1  g  of  initial  solid  explosive 

0^  specific  heat  at  constant  value. 

Subscripts: 

8  =  solid;  g  =  gas;  S  =  isentroplc;  H  =  dynamic  adiabatic 
(Hugonlot);  1  =  isothermal;  o  -  initial  value. 
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Chapter  5 

MATHEMATICAL  MODEL  PROPOSED 


5.1.  PRESENTATION  0?  THE  MODEL 

In  order  to  assure  ourselves  of  a  better  representation 
of  the  experimental  phenomena,  we  Btudied  a  mathematical  model 
which  takes  Into  account  the  remark's  made  above.  In  addition 
we  have  tried  to  re-obtain,  by  calculation,  some  experimental 
results  wnlch  are  not  obtained  by  the  current  models. 

5.1.1.  System  of  Equations 

The  solid  explosive  and  the  detonation  products  are 
treated  as  homogeneous  bodies,  and  it  13  assumed  that  the 
reaction  goes  to  completion: 

Solid  explosive  — »  detonation  products 

We  have  chosen  two  separate  equations  of  state  for 
representing  the  solid  and  the  gas,  but  we  assumed  that  at 
every  moment,  In  a  given  mesh,  they  are  at  the  same  pressure 
even  though  their  temperatures  are  different.  These  hypotheses 
have  been  discussed  In  the  preceding  chapter.  (For  notations, 
see  table  at  the  end  of  the  preceding  chapter.) 

For  the  solid,  we  have  tak>n  an  equation  of  state 
having  the  form: 


(1) 


where  Is  the  isothermal  pressure  (a  function  of  *  alone) 
which  we  have  chosen  as  given  by  the  equation  of  Pack,  Evans 


P* 


g(T  -  T0)  p,  *  j  exp  [p(,  .  ?)>  J 
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and  James,  and  g  is  the  coefficient  <**)  ,  which  Is  assumed  to 

v 

be  constant.  The  change  of  internal  energy  from  conditions 
p0,  vn,  T0  to  conditions  p,  v,  T  is  then  furnished  by 

e.  -  E..  «  C%  (T.  -  T„>  4|i;,  (v,  .  v„)  -  jr*.  pt  dv.  (2) 


Pbr  the  gas,  Cook's  equation  [4,  5]  was  employed,  even 
though  a  priori  it  had  been  determined  for  a  quite  different 
application.  In  effect,  this  law  tries  to  express  the  state 
of  the  detonation  products  under  the  conditions  of  Chapman- 
Jouguet  ( C . J .  )  with  the  reaction  terminated,  independently  of 
their  composition,  while  here  we  are  looking  for  a  law  which 
makes  it  possible,  for  a  given  explosive,  to  follow  the 
behavior  of  the  detonation  products  during  the  reaction. 
Nevertheless,  since  the  products  formed  are  essentially  the 
same  regardless  of  the  explosive  being  investigated,  we  are 
Justified  in  believing  that  the  error  made  by  choosing  Cook's 
law  is  quite  small. 

p  [v  -  ot (v) ]  •  n  R  T  (3) 

The  use  of  an  equation  of  state  of  the  form  p  =  RT  f(v) 
leads,  for  the  determination  of  the  variation  of  the  Internal 

energy,  to  the  simple  equation  (Tlx  •  p  m  °) 


(4) 


The  specific  heat  of  the  gas  varies  considerably  with 
the  temperature.  It  may  be  estimated  that  it  goes  from  0.1 5 
cal/g/°K  at  300^  to  0.8  cal/g/°K  at  4,000r)X.  Hence  an 
average  value  of  0.4  cal/g /°X  may  be  accepted,  all  the  more  so 
since  It  leads  to  quite  satisfactory  values  for  the  C.J. 
characteristics  of  the  explosive. 

The  knowledge  of  the  Internal  energy  of  the  solid  and 
the  gas  then  makes  it  possible  to  determine  that  of  the  mixture. 
By  calling  m  the  mass  of  gas  formed  at  a  given  instant  per  gram 
of  Initial  solid  explosive,  this  energy  Is  given  by 

E  -  E.  -  m(E,  -  E„)  Ml  -  m)  (E,  -  E„)  (5) 
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by 


Likewise,  the  specific  volume  of  the  mixture  is  defined 


v  *  m  v,  ♦  (1  -  m)v, 


(6) 


The  equation  of  chemical  kinetics  giving  the  rate  of 
liberation  of  chemical  energy  was  chosen  in  various  ways  so  as 
to  show  the  difference  between  Arrhenius'  law  and  the  linear 
law  as  a  function  of  pressure.  The  latter  has  even  been  some¬ 
what  complicated  so  as  to  account  for  certain  experimental 
results:  its  formulation  will  be  specified  during  the  discussion 
of  the  theoretical  results  obtained  by  this  mathematical  model; 
in  a  general  fashion  we  shall  write  it  as: 


f(p.  t.  ...) 


(7) 


Finally,  if  the  viscosity  and  heat  conduction  are  neg¬ 
lected,  the  fundamental  hydrodynamic  equations  furnish  the 
following  equations  expressed  in  Lagrange  variables  (X,  t): 


22L.J l 

dx  V 

(conservation  of  mass) 

(8) 

D  i£  .  o 

p*  3t  3X 

(momentum ) 

(9) 

is. .  Q  42  +  pi>X.  0 

at  Wr  a  t  p  at 

(energy  equation) 

(to) 

ax 

— —  »  u 

at 

(material  velocity) 

(ID 

In  this  way  we  have  a  system  of  eleven  equations  for 
the  twelve  functions  of  the  two  variables  X  and  t: 

u,  x.  v,  m,  E,  p,  V,,  T,  ,  E, ,  vf.  Tf,  Ef 

Hence  we  have  the  choice  of  an  additional  equation. 

We  shall  assume,  since  the  solid  cannot  heat  up  by  thermal  , 
conduction,  that  its  behavior  is  linked  solely  to  the  pressure 
variations  of  the  mixture.  Hence  It  can  only  be  subjected  to 
lsentroplc  transformations.  Since  we  have  already  assumed 
that  the  solid  and  gas  are  at  the  same  pressure  at  all  times  -- 
in  a  grid  --  (pB  =  pff  =  p),  the  equation  sought  is  therefore 

d  E,  +  p  dv,  •  0  (12) 
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For  the  sake  of  facility  of  use  in  the  subsequent  calculations 
we  have  replaced  it  by 


T. 


•  T>  tsp 


<v.. 


(12) 


This  system  of  twelve  equations  is  obviously  used  only 
during  the  reaction,  when  solid  and  gas  are  both  present.  It 
reduces  to  a  system  of  six  equations  in  the  following  two 
cases : 


The  reaction  is  not  initiated,  m  =  0,  the  functions 
are  then:  u ,  x ,  v  =  va ,  E  =  Es ,  p,  Ta  and  the  system  is  made 

up  of  equations  1,  2,  8,  9,  10,  11. 

The  reaction  ie  terminated,  m  =  1 ,  the  functions  are 
then:  u,  x,  v  =  v£,  E  =  E?,  p,  T  .  The  system  is  then  made  up 

of  equations  3,  A,  8,  9,  10,  11. 

The  impact  is  furnished  by  a  metal  plate  of  thickness  e, 
moving  with  a  velocity  V^.  The  pressure  at  all  points  of  the 

plate  is  assumed  to  be  zero. 

5.1.2.  Conduct  of  the  Numerical  Calculation 

The  choice  of  the  system  of  equations  of  finite  differ¬ 
ences  1 s  directly  Inspired  by  that  of  Rlchtmyer  [12]  with 
pseudoviscoslty  (q). 

The  three  hydrodynamic  equations  8,  9,  11,  are  then 
written  as: 


v..i 


-  i 

J± 


«.> 

1± 


-  (p  ♦  q)’  +  (p  ♦  ql'J.! 
P,  AX 


(8') 

(9') 

(11’) 


The  pseudoviscosity  is  given  by  the  equations 
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where  1  is  a  length. 

The  knowledge  at  time  n  of  the  functions  p,  v,  q  in  the 
grids  J  -  1,  J,  J  ♦  1  and  of  functions  u,  x  at  the  interfaces 

J  -  J  ♦  makes  it  possible  to  calculate  at  time  n  ♦  1  by 

means  of  equation  9':  by  means  of  equation  11*:  and 

by  means  of  equation  8':  vj“. 

Starting  from  m^,  equation  7  furnishes  m”4^,  regardless 
of  the  law  of  chemical  kinetics  employed. 


with 


*»;“  •  ntj  ♦ 
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rip.  T,  ...) 


(7’) 


Thus  we  have  eight  equations  left  for  the  determination 
of  the  eight  functions  E,  p,  v8 ,  T0 ,  E0 ,  v  ,  T  ,  E».  The  solu- 

tlon  of  the  system  of  eight  equations  with  eight  unknowns  is 
carried  out  in  the  following  manner:  equations  5  and  10  are 
calculated  as  a  function  only  of  the  variables  v8  and  v^,  them¬ 
selves  related  by  equation  6.  Hence  it  is  possible  to  determine 
by  iteration  the  values  of  these  variables  which  simultaneously 
satisfy  the  three  equations  5»  6  and  10.  Hence  this  calculation 
makes  it  possible  to  determine,  in  a  stepwise  manner,  the 
behavior  during  the  reaction  of  the  mixture  of  explosive  solid 
and  detonation  products  in  space  and  t.me.  It  should  never¬ 
theless  be  noted  that  equation  10  involves  the  pseudoviscosity, 
but  this  does  not  essentially  complicate  the  solution  of  the 
system  of  equations. 
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In  all  these  equations  subscript  J  has  been  omitted. 

During  the  moments  when  only  the  solid  exists  (m  =  0) 

In  a  grid,  or  the  solid  Is  completely  converted  Into  gas 
(m  =  1),  the  mathematical  model  Is  reduced  to  six  equations. 
Three  of  them,  which  are  fundamental  hydrodynamic  equations, 
Immediately  furnish  u,  x,  and  v  (va  or  vg)  at  moment  n  ♦  1 . 

Then  the  system  Is  reduced  to  three  equations  comprising  the 
three  unknowns  p,  E_  (or  E  ) ,  T.  (or  T^). 

s  go  D 

The  Initiation  of  the  reaction  was  carried  out  as  fol 
lows:  During  the  rise  of  the  shock  front  --  a  gradual  rise  c 
to  the  use  of  pseudoviscosity  —  the  explosive  Is  treated  as 
an  Inert  substance,  and  the  reaction  Is  released  only  when  the 
pressure  maximum  due  to  the  Impact  has  been  attained.  This 
state,  situated  on  the  Hugonlot  of  the  Inert  solid,  denoted  by 
subscript  H,  constitutes  the  Initial  conditions  for  the  deter¬ 
mination  of  the  subsequent  behavior  of  the  solid  during  the 
reaction  (lsentroplc  transformation).  Moreover,  the  tempera¬ 
ture  Tu  Is  taken  as  an  Initial  value  In  the  two  equations  of 
the  determination  of  the  Internal  energy  of  the  gas,  T  °  =  T„ 

6  M 

(equation  A),  and  of  the  chemical  kinetics  in  the  case  where 
the  latter  is  Arrhenius'  equation. 

5.1. 3.  Choice  of  Numerical  Values 

The  object  of  this  study  Is  more  that  of  determining 
the  Influence  of  the  various  parameters  on  the  generation  of 
the  detonation  In  an  explosive  than  to  explain  the  particular 
quantitative  behavior  of  one  of  them.  Nevertheless,  we  have 
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made  an  attempt  to  choose  the  various  constants  In  such  a  way 
as  to  represent  explosive  D  in  a  satisfactory  manner  [2]. 

For  the  equation  of  state  of  the  solid  explosive,  the 
experimental  measurements  of  the  dynamic  adiabatic,  carried 
out  on  this  explosive  l8]  gave  the  following  linear  law  U(u): 

u  *  2  400  ♦  i,66  u  (mts)  ( meter-  ton- second ) 

The  value  of  the  specific  heat  of  this  explosive  [ll] 
has  led  us  to  choose  Cy,  which  we  assume  constant  regardless 

of  the  temperature,  as  equal  to  880  kj/ton. 

On  the  basis  of  these  results  the  constants  of  the  PEJ 
equation  were  adjusted,  and  the  values  retained  are  as  follows: 


a  *  1,349.10*  p  ■  16.5  g  .  3.491  (meter- ton-second ) 
for  an  initial  density  and  temperature  of 

p.  •  1.70  t/ro>  i;  •  300*  K 


respectively. 

By  means  of  the  equations  given  in  Chapter  1 ,  it  is  now 
possible  to  determine  the  various  functions  necessary  for  the 
calculation  of  the  behavior  of  the  inert  explosive  during  the 
shock-detonation  transition,  in  particular  the  dynamic  adiabatic 
of  pole  p  =  0,  v  =  vn ,  and  the  lsentroplc  lines  each  originating 

from  a  point  of  this  adiabatic. 

The  determination  of  the  detonation  characteristics  of 
this  explosive  then  consists  of  solving  the  system  of  six 
equations  recapitulated  below: 


p  (v  .  tt)  *  n  r  t  (equation  of  state  ) 
u  .  d(i  - (continuity  equation) 

\  V#  / 


d’  « 


(momentum  equation) 


A  E  -  Q,  ■  |  p  (v.  -  v) 


(energy  equation) 
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4*  •  C,  (T  ^  TJ 
D  -  u  •  • 


(equation  of  Internal  energy) 
(Chapman- Jouguet  condition) 


The  solution  of  this  system  of  equations  may  be  carried 
out  in  the  following  manner  [2]: 

The  determination  of  the  velocity  of  sound  a^  =  (^)B 

along  an  isentroplc  dE  ♦  p  dv  =  0  with  dE  =  CydT,  talcing  into 
account  the  equation  of  state  selected: 


(v  -  a )  dp  ♦  p  ( l  -  dv  ■  n  R  d 


is  furnished  by 


&).  ■-  <■-&♦!*) 


v  -  a 


or 


0 


da  n  R\  p  v» 
dv  C,/  v  -  a 


Substituting  this  value  into  the  C.J.  condition,  after 
replacing  the  lefthand  side  of  the  latter  by  the  values  of  D 
and  u  taken  from  the  continuity-  and  momentum  equation,  we  get 

.  da  .  n  R  v  -  a 

1  -  “T“  ▼ 

dv  C,  v.  -  v 


The  use  of  the  law  r*.  (v) ,  given  by  v  •»,«'♦  b,  «♦  c,  , 

yields 


ia. _ 1 - 

dv  2  a  ♦  bj 


The  equation  to  be  solved  then  reduces  to 


l  ♦ 


n  R 
C, 


▼  -  a 
v,  -  v 


1 


The  lefthand  side  of  this  equation  is  constant,  hence 
v  and  ot are  obtained  directly.  Tne  combination  of  the  energy 
equation  leads,  after  elimination  of  the  temperature,  to  the 
determination  of  p  by  means  of 
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p  [(v  *  a>  *  7^  (v»  ‘  v|]  "  '  R  [T.  +  -§:-J 

The  other  equations  then  permit  solving  the  problem 
completely. 

V/e  have  chosen  the  function  ot(v),  directly  Inspired  by 
Cook's  curve  [4,  5 J »  given  by 

v  •  a,  625  a*  -  0,8  a  ♦  0,34 

Using  the  following  numerical  values  (YTS) 

n  •  32.108  R  ■  B.Suse.lO*5  C,  «  1.672  Q,  •  5.225. 1^ 

T,  ■  300*  K 


calculation  yields 

p  ■  216.10’  pitzes  D  *  8.305  rr> /«*  u  •  1.529  m/s«C 
T  »  4. 124*  K 

f^nly  the  velocity  of  detonation  Is  slightly  too  high 
with  regard  to  the  values  of  pressure  and  material  velocity. 

A  more  elaborate  equation  of  state  would  permit  a  better  adjust¬ 
ment  of  these  values  with  respect  to  each  other,  but  apart  from 
the  fact  that  It  would  considerably  complicate  the  calculation, 
It  would  probably  contribute  only  a  very  slight  Improvement  to 
the  proposed  scheme. 

5.2.  RESULTS 

This  system  of  equations  was  adapted  for  numerical  cal¬ 
culation  on  an  electronic  digital  computer.  Its  stability 
was  verified  numerically.* 

5.2.1.  Arrhenius'  Law 

A  first  series  of  calculation  was  carried  out  using 
Arrhenius'  Law  as  the  law  of  chemical  kinetics: 

dm  . .  , 

—  ■  v  (i  -  m)  e‘«7 


♦The  program  was  prepared  by  Yr.  Vldart  and  Yme  Frouteau, 
and  used  on  the  Bull  Gamma  60  electronic  computer. 
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The  temperature  T  was  chosen  as  thet  of  the  gas  (Tg). 
However,  the  use  of  two  different  temperatures  (Ta  and  Tg)  for 

the  same  pressure  (the  gas  is  at  A.OOO^K  when  the  solid  is  at 
1,800°K)  entails  the  rapid  initiation  of  the  reaction.  Then 
the  latter  is  complete  in  a  short  time,  since  Arrhenius'  law 
is  very  sensitive  to  the  temperature. 

Under  these  conditions  the  detonation  starts  at  the 
Interface  as  in  the  case  of  homogeneous  explosives  but  takes 
place  so  rapidly  that  there  is  hardly  any  "incubation"  period, 
since  the  inert  shock  does  not  have  the  time  to  be  propagated 
in  the  explosive. 

Nevertheless,  there  is  obtained  a  gradual  rise  of  the 
velocity  of  the  shock  front  which  3,8  reactive  from  the  begin¬ 
ning  (the  pseudoviscosity  is  only  an  artifice  for  the  calcu¬ 
lation),  but  the  reaction  is  such  that  the  steady  detonation 
state  is  attained  very  fast  and  without  the  possibility  of  a 
slowdown. 

5.2.2.  The  Law  of  Pressure 

The  following  numerical  study  was  not  oriented  toward 
the  explanation  of  the  quantitative  behavior  of  a  particular 
explosive  but  rather  toward  determining  the  influence  of 
various  parameters  on  the  process  of  initiation  of  the  detona¬ 
tion,  in  order  to  compare  their  variation  with  the  experimental 
values . 


In  effect,  the  use  of  a  numerical  method  with  finite 
differences  for  the  solution  of  the  system  of  equations  is 
not  without  introducing  some  inconvenience.  In  order  to 
reduce  the  fluctuations  as  much  as  possible,  we  have  been  led 
to  choosing  a  dense  grid  in  space  and  time.  The  grid  chosen 
(/lx  =  0.1  mm ,  At  -  5  n  sec)  requires  a  large  number  of  calcula¬ 
tions.  So  as  to  make  sure  that  a  complete  calculation  does  not 
require  a  prohibitive  number  of  machine  hours,  we  were  forced  — 
while  preserving  an  order  of  magnitude  compatible  with  the 
experimental  values  --  to  obtain  the  stable  detonation  state 
in  approx.  2  jusec  for  the  "averaee"  case  (Impact  velocity  of 
a  copper  projectile:  1,600  m/sec),  even  though  this  is  slightly 
low  in  absolute  value. 

Since  we  wich  to  study  the  behavior  of  granular  ex¬ 
plosives,  we  were  first  motivated  by  Eyrlng's  law  grain-burning 
L 7 3 »  then  by  the  law  of  holeburning. 


93 


In  the  first  law  it  la  assumed  that  the  exploalve  ia 
preaent  in  the  form  of  apherea  having  an  average  radius  R. 

The  chemical  decomposition  takes  place  only  on  the  surface, 
and  the  Influence  of  the  pressure  on  the  rest  of  the  grain  ia 
neglected.  Strictly  speaking,  with  the  characteristics  which 
we  have  chosen  for  the  solid  explosive,  when  the  pressure 
changes  from  100  kb  to  200  kb  the  radius  of  the  grain  decrease 
by  about  3.5,2. 

This  hypothesis  leads  us  to  the  assumption  that  the 
mass  of  the  sphere  of  radius  R  is  directly  proportional  to  its 
volume  during  the  entire  reaction. 

By  calling  m  the  mass  of  gas  formed  per  gram  of  initial 
explosive,  and  r  the  radius  of  the  sphere  at  moment  t, 

r>  , 

F  '  1  '  m 


whence  we  get:  r  *  R  u  . 
and  by  differentiation: 


3  r'  dr  p  dm 
R>  d  t  dt 


Hence  the  equation  of  the  rate  of  reaction  has  the 
following  form: 


dm 

dt 


In  this  case,  dr/dt  is  negative,  and  the  radius  of  the 
grain  decreases  during  the  reaction. 

,  2 

The  terra  (1  -  m)j  is  the  equivalent  of  the  "geometric 

form  function"  of  the  combustion  equations  of  powders,  while 
dr/dt  corresponds  to  the  linear  combustion  rate  and  dm/dt  to 
the  rate  of  consumption. 

Next,  our  reasoning  was  influenced  --  under  the  same 
assumptions  —  by  holeburning.  Then  the  reaction  is  assumed 
to  propagate  in  a  divergent  spherical  manner  from  the  center 
of  the  explosive  grain.  We  then  get 


j±  t 

R> 


m 
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whence 

i*  ■  R  mf 

and 

8  r*  _d”  dm 

RJ  d»"  at 

whence 

jfel  .  _L  dr 

d»  R  *  *  ’  *  dt 

this  time  dr/dt  being  positive. 


CoabuttiM  Mr  fraln  (J)  CortMtlon  Mr  troa  <$ 

(arain-b«rnlna)  (ho!#Mrnlna) 

(1)  -  direction  of  propagation  of  the  reaction;  (2)  - 
grain-burning;  (3)  -  holeburning. 

The  essential  difference  between  these  two  equations  of 
chemical  kinetics,  all  things  being  equal,  is  that  the  first 
leads  to  a  slowing  down  of  the  rate  of  reaction  dm/dt  in  direct 
proportion  to  the  increase  in  the  mass  (m)  of  gas  formed,  while 
the  second,  on  the  contrary,  causes  its  acceleration. 

The  last  point  which  we  have  to  make  explicit  for  carry¬ 
ing  out  the  complete  calculation  is  the  manner  in  which  the 
radius  of  the  grain  varies  in  time,  in  other  words,  specify 
the  combustion  rate  dr/dt. 

ThlB  variation  depends  on  the  thermodynamic  conditions 
of  the  grain  and  of  the  medium  surrounding  it.  It  is  not 
possible  at  the  present  time  to  determine  it  experimentally, 
as  a  function  of  the-  various  parameters  (pressure,  temperature, 
etc.).  We  can  only  be  guided  in  our  choice  by  analogies  with 
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similar  phenomena  (chemical  reactions,  combustion)  and  the 
comparison  between  the  theoretical  and  experimental  results: 

--  The  predominant  influence  of  pressure  in  the  tran¬ 
sition:  deflagration-detonation  (see  Chapter  2); 

—  The  progress  diagram  of  the  reactive  shock  front 
during  the  initiation  of  the  detonation  In  a  granular  explosive 
where  the  compression  waves  accelerate  the  shock  as  they  over¬ 
take  It; 


—  The  combustion  of  powders  both  at  low  pressure  (solid 
propergols )  and  at  several  thousand  atmospheres  (lnterml 
ballistics  of  cannons); 

--  The  measurements  made  by  certain  researchers  up  to 
10  kb  on  solid  explosives  L 9 ] . 


The  whole  body  of  these  partial  data  has  led  us  to  prefer, 
as  Independent  variable,  the  pressure  to  the  temperature,  all 
the  more  so  since  It  is  directly  involved  In  the  hydrodynamic 
equations. 

Hence  we  are  led  to  choose  a  law  of  the  form:  |-^|  .  ap«  +  b 
with  a  close  to  1 . 

In  first  approximation,  we  chose 


a 


1 


or 


I- 1 

'dt  1 


•p  ♦  b 
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The  use  of  a  linear  law  le  a  priori  quite  satisfactory, 
since  we  have  seen  that  the  choice  of  an  exponential  law  entails 
a  much  too  abrupt  release  of  the  reaction,  and  does  not  allow 
the  representation  of  the  experimental  results  of  solid  granular 
explosives . 

The  majority  of  theoretical  results  presented  below  were 
obtained  by  means  of  this  linear  law  of  pressure. 

Under  the  assumption  of  a  reaction  governed  by  the  law 
Inspired  by  grain-burning,  the  equation  of  chemical  kinetics 
defined  above  Is  written  as 

“jjp-  (i  *  Up  +  b) 


In  order  to  Ta^e  evident  the  pressure  threshold  p^  below 

which  It  seem9  that  the  reaction  cannot  be  Initiated  experi¬ 
mentally,  we  gave  b  the  value  of  (-  ap^). 

Then  the  final  equation  employed  Is 


m) 


(p  -  p^ 


and  under  the  assumption  of  the  law  Inspired  by  holebuming,  we 
get 


-  Diagram  of  the  progress  of  the  reactive  shock  front: 

The  use  of  a  law  of  chemical  kinetics  depending  on  the 
pressure,  combined  with  a  function  having  a  form  analogous  to 
that  of  grain-burning,  furnishes  a  progress  diagram  of  the 
reactive  shock  which  shows  a  rather  good  agreement  with  experi¬ 
mental  results.  In  effect  we  obtain,  Immediately  upon  the 
passage  of  the  shock,  a  very  high  rate  of  reaction.  Hence  the 
pressure  Increases  abruptly,  end  then  gradually  decreases  to 
the  conditions  of  Chapman- Jouguet.  The  reactive  s..ock  front 
accelerates  from  fhe  moment  of  Impact,  to  attain  gradually  the 
steady  detonation  state  (diagrams  1  and  J>).  By  contrast,  In 
the  case  of  holeburning,  the  reaction  Is  quite  weak  after  the 
passage  of  the  shock.  The  latter  is  therefore  propagated  In 
the  explosive  at  an  almost  constant  rate,  as  In  an  Inert 
medium.  It  Is  at  the  target-projectile  interface,  the  region 
which  le  compressed  for  the  longest  period  of  time,  that  the 
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reaction  first  becomes  considerable;  then  compression  waves 
are  propagated  toward  the  shock  front,  the  pressure  signal 
gradually  straightens  out  and  finally  assumes  a  form  analogous 
to  that  obtained  under  the  assumption  of  grain-burning.  Hence 
during  the  first  moments  the  reactive  shock  front  propagates 
at  a  constant  velocity,  then  it  becomes  accelerated  quite 
abruptly  when  the  compression  waves  originating  from  the  inter¬ 
face  catch  up  with  it,  and  finally  gradually  attains  the 
velocity  of  the  steady  detonation  state.  Hence  a  slight  slope 
discontinuity  is  observed  in  this  case  in  the  progress  diagram 
(x,  t)  of  the  shock  in  the  explosive  (diagrams  2  and  3)* 

The  use  of  a  form  function  involving  the  average  diameter 
of  the  grains  makes  it  possible  to  represent,  at  least  quali¬ 
tatively,  the  following  experimentally  observed  finding:  the 
smaller  tne  average  grain  diameter,  the  easier  it  is  to  bring 
about  the  detonation,  and  the  faster  the  steady  state  is 
attained.  It  should  nevertheless  be  pointed  out  that  the 
threshold  of  the  initiation  of  the  reaction  does  not  intervene 
in  this  model  simultaneously  with  the  radius  of  the  grain,  as 
may  be  assumed  on  the  basis  of  experiments. 

Three  calculations  were  carried  oui  for  grain  sizes  of 
50,  ICO  and  140  for  the  case  of  grain-burning.  The  results 
are  shown  on  diagrams  1  (R  =  140  y),  4,  5  and  6. 

If  we  consider  the  experimental  results  in  greater 
detail,  the  progress  diagram  of  the  shock  front  exhibits  two 
Important  characteristics: 

--  The  acceleration  takes  place  from  the  moment  of 

impact ; 


--  There  is  an  abrupt  change  in  this  acceleration. 

The  above  calculations  do  not  lead  to  these  results. 

The  law  of  reaction  inspired  by  grain-burning  expresses  only 
the  first  of  these  characteristics,  while  the  law  inspired  by 
holeburning  expresses  only  the  eecond. 

To  obtain  the  acceleration  of  the  shock  front  from  the 
m<  „ent  of  Impact,  It  is  necessary  that  the  reaction  br  initiated 
from  this  moment  to  a  non-negliglble  extent.  This  is  realized 
by  the  uee  of  a  gradual  law  of  energy  liberation,  ior  example, 
a  linear  law  as  a  function  of  the  pressure,  which  is  more 
gradual  than  any  exponential  law. 
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The  abrupt  change  of  acceleration  can  be  conceived  only 
if  an  abrupt  excess  pressure  appears  behind  the  shock  front. 
This  excess  pressure  can  be  produced  only  by  the  arrival  of  a 
train  of  compression  waves  (or  a  second  shock)  which,  arising 
in  the  already  perturbed  explosive  zone  between  the  impact 
surface  and  the  front  of  the  first  shock,  overtakes  the  latter. 

This  is  the  explanation  given  for  interpreting  the 
experimental  results  obtained  with  the  homogeneous  explosives 
where  it  is  assumed  that  the  second  shock  arises  at  the 
barrier-explosive  Interface  (see  preceding  chapter). 

However,  in  the  case  of  the  granular  explosives,  the 
second  shock  ought  to  be  propagated  in  a  medium  which  has 
already  strongly  reacted  --  the  progress  diagram  (x,  t)  of  the 
Initial  shock  3hows  an  acceleration  from  the  moment  of  impact. 
Moreover,  the  reaction,  after  the  passage  of  this  wave  train, 
does  not  have  to  be  complete,  since  the  velocity  of  the  steady 
detonation  state  is  attained  only  gradually,  and  by  a  lower 
value  after  the  slope  discontinuity  of  the  progress  diagram. 

In  order  not  to  obtain  excess  velocity  --  as  in  the 
case  of  the  homogeneous  explosives  --  before  the  establishment 
of  the  steady  detonation  state,  we  have  been  led  to  believe 
that  it  is  in  the  vicinity  of  the  initial  shock  front  that  the 
excess  pressure  is  brought  about  abruptly.  Accordingly  this 
would  be  provoked  by  an  abrupt  liberation  of  chemical  energy, 
and  hence  would  correspond  to  a  Jump  in  the  evolution  of  the 
"reaction."  This  abrupt  liberation  of  energy  may  be  accepted 
if  we  remember  that  several  reactions  are  involved,  and  that 
the  chemical  kinetics  of  the  whole  phenomenon  is  not  known. 

This  assumption  is  strengthened  by  the  following  obser¬ 
vations: 


--  The  examination  of  the  decompo c ltion  of  metal 
nitrides  has  revealed  an  abrupt  variation  of  the  reaction  rate 
as  a  function  of  temperature  t l] . 

--  The  study  of  electric  conductivity  in  the  reaction 
zone  of  various  granular  explosives  (detonation  in  steady 
state)  reveals  notable  differences  in  the  form  of  the  electric 
signal  recorded.  The  observations  make  one  assume  that  the 
appearance  of  ionization  is  Intimately  linked  with  ‘.he  chemical 
reaction.  The  deviations  seem  to  be  due  to  the  differences 
with  respect  to  the  mechanism  of  reaction,  with  the  transitory 
appearance  of  Intermediary  products,  related  to  the  physical- 
chemical  properties  of  the  explosive  in  question,  while  the 
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determination  of  the  detonation  characteristics  with  the  aid 
of  the  hydrodynamic  theory,  where  it  Is  assumed  that  a  thermo¬ 
dynamic  equilibrium  has  been  achieved,  only  considers  the  final 
products  L 1 3] - 

--  The  behavior  of  certain  explosives  where  a  state  of 
"low  order  detonation"  is  established,  gives  rise  to  the  belief 
that  the  liberation  of  chemical  energy  is  capable  of  not  being 
total  under  certain  conditions,  while  still  assuring  a  steady 
state  tlA], 


d'lprAl  AUDUBLRT  -  Rdf  1  l£)  ©  d'aprt,  R.  SC MALL  •»  *.  VOLLRATH.  R4f  13 


(1)  -  rate  of  reaction  (arbitrary  units);  (2)  - 
according  to  Audubert  [ 1 ] »  (3)  -  conductivity; 

(A)  -  according  to  R.  Schall  and  K.  Volirath  [13]» 

(5)  -  probe. 

It  follows  from  these  remarks  that  the  mechanism  of  the 
"reaction"  depends  on  the  explosive  In  question  and  varies 
during  the  evolution  all  the  chemical  reactions  which  lead 
to  detonation. 

This  explains  why  It  Is  sometimes  possible  to  observe, 
experimentally,  an  "overshoot"  during  the  establishment  of  the 
steady  state  of  detonation  of  certain  explosives  lIO].  The 
ease  with  which  the  train  of  compression  waves  Is  transformed 
into  a  shock  is  greater  when  the  Intensity  of  the  initial 
shock  Is  lower  —  a  longer  time  for  the  establishment  of  this 
shock  by  the  combination  of  the  various  compression  waves  which 
arise  In  a  medium  which  has  not  reacted  to  a  great  extent. 

Hence  the  overshoot  Is  greater  when  the  Initial  shock  Is 
weaker  (see  Jacobs  diagram.  Chapter  A). 
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Hence  the  general  case  —  without  overshoot  —  requires 
the  introduction  of  a  discontinuity  into  the  liberation  of 
chemical  energy  in  the  vicinity  of  the  shock  front. 

The  use  of  an  exponential  law  gives  a  good  overall 
representation  of  a  threshold  state,  but  we  have  already  seen 
that  the  “discontinuity"  Introduced  in  this  fashion  is  much 
too  abrupt  to  represent  the  real  situation. 

The  use  in  the  calculation,  presented  above,  of  a  form 
function  analogous  to  holeburning  does  entail  a  progressive 
acceleration  of  the  reaction  in  proportion  to  its  state  of 
advance,  but  the  slope  discontinuity  obtained  in  this  way  is 
small,  while  the  shock  propagates  during  the  first  moments  as 
an  inert  shock. 

While  a  combination  of  the  two  laws  of  reaction,  having 
the  form 


dm 

dt 


3  a  , 

~ r-  fci  n*  ♦  k, 
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indeed  furnishes  both  the  acceleration  from  the  moment  of 
impact  and  the  velocity  Jump  obtained  experimentally,  these 
two  characteristics  are  much  too  attenuated  to  permit  us  to 
consider  this  law  of  chemical  kinetics  satisfactory.  In  effect, 
each  mode  of  reaction  weakens  the  peculiarity  contributed  by 
the  other.  In  this  equation  coefficients  kj  and  k2  would 

represent  the  percentages  of  each  of  the  two  geometric  forms, 
values  linked  to  the  charge  density  of  the  explosive. 

Calculation  leads  to  the  following  result:  The  pressure 
gradually  rises,  in  the  course  of  time,  in  the  zone  perturbed 
by  the  shock,  while  at  a  given  Instant  the  preesure  maximum 
i 8  situated  either  at  the  target/pro Jectlle  interface  or 
immediately  behind  the  shock  front,  depending  on  the  relative 
values  of  kj  and  k2.  In  all  cases  the  deviation  is  email,  and 

the  progress  diagram  of  the  shock  front  is  intermediate  between 
those  obtained  above  for  each  of  the  two  modes  of  reaction  (by 
grain  or  by  hole). 

In  view  of  the  lack  of  success  with  these  different 
methods,  we  have  chosen  the  model  of  grain-burning  --  which 
furnishes  the  acceleration  of  the  shocn  front  from  the  moment 
of  impact  —  and  the  slope  discontinuity  in  the  progress 
diagram  of  the  shock  was  obtained  by  the  artificial  introduc¬ 
tion  of  a  Jump  in  the  liberation  of  the  chemical  energy.  This 
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Jump  was  related  to  a  threshold  value  of  the  pressure  since 
at  a  given  moment,  the  pressure  Is  maximum  In  the  vicinity 
of  the  shock  front  in  this  model  (Diagram  7). 

-  Threshold  of  Initiation  of  Detonation: 

Since  the  pressure  threshold  is  made  directly  explicit 
in  the  equation  of  chemical  kinetics,  it  is  a  known  quantity. 
Hence  the  scueme  presented  here  does  not  permit  ltB  theoretical 
determination. 

-  Initial  Point  of  Release  of  the  Detonation: 

The  experimental  results  show  that  the  detonation  is 
initially  established,  In  the  case  of  granular  explosives,  in 
the  interior  of  the  explosive  and  not  on  the  impact  surface 
as  is  the  case  for  the  homogeneous  explosives. 

The  use  of  a  very  gradual  reaction  velocity  makes  it 
possible  to  bring  about  a  sufficiently  slow  rise  in  the  reac¬ 
tion  rate  so  that  the  hydrodynamic  equations  may  intervene  in 
the  course  of  the  reaction.  The  gradual  pressure  rise  on  the 
impact  surface  of  the  explosive  is  slowed  down  by  the  pertur¬ 
bations  produced  by  the  equality  of  pressures  which  is  estab¬ 
lished  in  this  plane  at  all  moments.  These  perturbations  are 
constituted  of  a  train  of  compression  waves  which  propagate 
in  the  projectile,  while  the  expansion  waves  propagate  in  the 
explosive . 

The  use  of  a  law  of  chemical  kinetics  linked  to  the 
pressure  then  makes  it  possible  to  observe  that  the  initial 
point  of  the  complete  reaction  (m  =  1  )  is  not  situated  at 
the  pro Jectlle/target  Interface  but  In  the  interior  of  the 
explosive. 

This  result  is  Btlll  sharply  different  from  that 
yielded  by  Arrhenius'  law  where  the  first  point  of  complete 
reaction  is  always  situated  on  the  impact  surface  of  the 
explosive . 

The  explanation  offered  by  various  researchers  of  this 
phenomenon  obtained  by  means  of  a  frontal  shock  (barrier- 
receiver  assembly)  as  being  the  consequence  of  expansion 
effects  due  to  the  edges  of  the  cartridge  of  explosive  is  in 
agreement  with  the  preceding  results.  Moreover  this  expansion 
is  not  t ne  only  one  which  Intervenes,  but  we  have  to  add  to  it 
that  which  rises  from  the  barrier.  Hence  these  two  results 
bring  about  a  decrease  of  pressure  on  the  surface  of  entry  of 


102 


the  explosive,  hence  a  decrease  of  the  reaction  rate,  since  we 
have  linked  the  latter  to  the  pressure.  Hence  the  mathematical 
model  proposed  here  makes  It  possible  to  explain  this  phenomenon. 

A  difficulty  nevertheless  persists,  since,  despite  the 
use  of  a  very  gradual  law  of  chemical  kinetics,  the  distance 
separating  the  Initial  point  of  detonation  (m  =  1 )  from  the 
surface  of  entry  of  the  shock  Into  the  explosive  is  quite 
small,  while  experimentally  the  distance  found  Is  much  greater. 

To  subject  this  problem  to  a  strict  treatment  It  would 
be  necessary  to  take  Into  account  the  two  expansions  Indicated 
above.  Then  the  phenomenon  becomes  two-dimensional  non- 
statlonary.  We  did  not  undertake  this  complex  study  with 
three  variables. 

However,  the  study  of  the  progress  diagram  of  the 
reactive  shock  front  has  shown  us  that  to  obtain  a  good 
agreement  between  theory  and  experiment  It  Is  necessary  to 
Introduce  a  discontinuity  In  the  chemical  kinetics.  This 
considerably  facilitates  then  the  obtainment  of  the  Initial 
point  of  detonation  In  the  Interior  of  the  explosive,  at  a 
considerable  distance  from  the  Impact  surface  (diagrams  7  and 
6). 


From  the  Initial  point  of  complete  reaction  the  end 
reaction  front  propagates  toward  the  shock  front  and  then 
reduces  the  thickness  of  the  zone  In  the  process  of  reaction 
until  the  establishment  of  the  steady  state. 

If  the  explosive  between  the  face  of  entry  of  the 
shock  and  the  Initial  detonation  front  has  reacted  only  slightly, 
which  favors  the  use  of  a  discontinuity  In  the  liberation  of 
chemical  energy,  the  phenomenon  of  "retonatlon"  Is  re¬ 
encountered  on  the  theoretical  level:  an  end  reaction  front 
propagates  toward  the  face  of  entry.  The  wave  associated  with 
this  phenomenon  Is  called  "retrograde"  wave. 

-  Influence  of  the  Intensity  of  tho  Pressure  Signal: 

The  intensity  of  the  pressure  signal  generated  In  the 
receptor  explosive  depends  on  the  material  of  the  projectile 
and  the  velocity  of  Impact.  Hence  the  modification  of  one  of 
these  two  parameters  makes  It  possible  to  vary  the  intensity 
of  the  pressure  signal  Induced  In  the  explosive  target.  We 
have  chosen  to  modify  the  Impact  velocity  by  using  a  copper 
projectile. 
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Three  calculations,  giving  respectively  65,  100  and  135 
kb,  have  been  carried  out.  They  made  it  possible  to  verify 
(diagram  9)  that  the  distance  covered  by  the  reactive  shock 
without  the  establishment  of  the  steady  detonation  state  is 
greater  when  the  pressure  is  lower.  This  steady  state,  attained 
in  an  asymptotic  manner,  does  not  permit  a  precise  determination 
of  the  "Induction  period,"  but  by  contrast,  the  velocity  jump 
in  the  progress  diagram  of  the  reactive  shock  is  easier  to 
demonstrate. 

These  results  are  in  rather  good  agreement  with  the 
experimental  results  [3]  since  they  furnish  a  gradual  variation 
of  thl9  "induction  period,"  in  contrast  to  Arrhenius'  law 
which  leads  to  a  pressure  threshold  below  which  detonation  does 
not  take  place,  and  above  which  the  time  of  establishment  of 
the  detonation  is  almost  independent  of  the  initial  conditions. 

A  calculation  was  attempted  for  an  induced  pressure  of 
180  kb:  numerical  instabilities  prevented  us  from  carrying  out 
this  study  successfully.  A  grid  more  compact  in  space  and  time 
for  the  solution  of  the  equations  with  finite  differences 
should  permit  solving  this  difficulty,  but  this  kind  of  modi¬ 
fication  was  not  undertaken  since  it  would  lead  to  very  lengthy 
periods  of  machine  calculation. 

-  Influence  of  the  Duration  of  the  Pressure  Signal 

The  experimental  studies  on  the  duration  of  the  pressure 
signal  as  a  function  of  the  thickness  of  the  projectile  have 
revealed  in  the  case  of  aluminum  that  the  decay  of  the  shock 
takes  place  much  sooner  in  space  and  time  than  could  be  pre¬ 
dicted  by  the  hydrodynamic  theory  [6].  To  explain  this 
phenomenon  it  is  suggested  that  the  elasto-plastlc  beha /lor  of 
the  material  constituting  the  projectile  cannot  be  neglected. 

To  our  knowledge,  analogous  experiments  have  not  yet  been 
carried  out  on  copper,  and  it  is  possible  that  the  hydrodynamic 
theory,  too,  i3  deficient. 

In  the  theoretical  study  which  we  have  undertaken  we 
did  not  take  into  account  this  characteristic,  and  we  have 
assumed  that  the  projectile  behaves  like  a  fluid.  On  the 
other  hand,  the  use  of  a  law  of  chemical  kinetics  which  is  a 
gradual  function  of  the  pressure  does  not  permit  stopping  the 
reaction  even  when  a  beam  of  expansion  waves  rises  in  the 
perturbed  products.  The  reaction  is  slower,  but  it  leads 
Inexorably  to  the  detonation  except  if  the  pressure  becomes 
less  than  20  kb,  the  threshold  which  we  have  chosen  for  its 
initiation . 
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Cn  Diagram  10  are  plotted  the  results  obtained  for  the 
two  projectiles  launched  at  the  same  velocity  (800  m/sec)  but 
having  different  thicknesses:  1  and  3  mm,  respectively.  At 
the  start  of  the  phenomenon,  the  two  curves  coincide  (0.5  fisec 
on  the  diagram),  since  the  beam  of  expansion  waves  which  rises 
in  the  projectile  has  not  yet  attained  the  interface.  They 
separate  at  around  0.6  jasec,  and  those  obtained  for  the  pro¬ 
jectile  of  1  mm  thickness  show  that  the  expansion  is  much 
greater  at  the  Interface  than  in  the  interior  of  the  explosive. 
At  2  jisec  the  detonation  develops,  leading  to  the  steady  state 
in  the  case  of  the  3  mm  projectile,  while  at  no  point  did  the 
reaction  rate  (m)  attain  the  value  of  1  in  the  case  of  the  1  mm 
projectile. 

-  Influence  of  the  Initial  Temperature 

To  determine  whether  the  proposed  model  is  sensitive  to 
the  initial  temperature  conditions  of  the  explosive,  two  cal¬ 
culations  were  carried  out.  Me  assumed  that  before  Impact, 
the  explosive  had  an  initial  temperature  of  300°K  in  one  case, 
350°K  in  the  other.  The  constants  of  the  equation  of  state 
of  the  inert  explosive  were  not  modified  for  such  a  small  tem¬ 
perature  variation  since  at  the  present  time  it  is  impossible 
to  distinguish  experimentally  between  the  two  dynamic  adlabatics 
corresponding  to  the  initial  conditions  of  300°K  and  350 nK. 

The  results  of  the  calculation  are  very  similar,  and 
exhibit  no  experimental  dispersion.  As  a  result,  in  this 
scheme,  the  temperature  is  only  of  slight  importance  which  is 
quite  different  from  the  result  obtained  by  Arrhenius'  law. 

Thus  the  agreement  with  the  experimental  results  is  quite 
satisfactory  L3J. 

-  Study  of  the  Steady  State  of  Detonation 

The  preceding  model  is  still  valid  when  the  stable 
state  of  detonation  has  been  est'^lisatso .  nence  it  furnishes 
a  schematic  representation  of  the  reaction  zone. 

In  the  first  place  It  should  be  noted  that  the  values 
of  the  detonation  characteristics  in  the  Chapman- Jouguet 
(C.J.)  plane  are  not  re- encountered  exactly  by  machine  calcu¬ 
lation.  In  particular,  the  detonation  velocity  is  considerably 
greater  (by  about  20*)  that  that  calculated  directly  by  means 
of  the  hydrodynamic  equations.  The  C.J.  point  itself  is 
slightly  displaced  in  the  ( p ,  v)  plane  with  respect  to  the 
theoretical  values.  It  is  possible  that  this  is  due  to  the 
introduction  of  the  pseudoviscosity  as  well  as  to  the  choice 
of  numerical  method  itself  (equations  with  finite  cif ferences ) . 
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The  examination.  In  plane  (p,  v),  of  the  behavior  of 
the  mixture  of  the  explosive  solid  and  detonation  products 
during  the  reaction  shows  (diagram  11)  that  it  starts  from 
point  H  of  the  Hugonlot  of  the  solid  defined  above,  first  In¬ 
creases  in  pressure,  then  decreases  to  attain  finally  the 
values  which  we  have  denoted  as  computed  C.J.  values,  rather 
close  to  the  theoretical  C.J.  conditions;  this  curve  is  inter¬ 
mediate  between  the  adiabatic  of  the  inert  solid  and  the 
isentropic  of  the  detonation  products  which  passes  through  the 
C.J.  point.  It  is  to  be  regretted  that  few  points  are  situated 
in  the  ascendant  part  of  this  curve  as  a  result  of  the  overly 
large  Ax,  At  grid  employed  in  the  equations  with  finite 
differences . 

It  may  be  observed  that  the  pressure  maximum  remains 
below  the  value  of  Von  Neumann's  peak.  This  is  normal,  since 
the  shock  is  immediately  reactive.  A  more  dense  grid  would 
perhaps  permit  to  reduce  even  further  the  pressure  maximum 
obtained  by  this  model. 

The  thickness  of  the  reaction  zone  depends  on  the  con¬ 
stants  introduced  into  the  law  of  chemical  kinetics.  Uith  the 
values  which  we  have  chosen,  It  is  of  the  order  to  one  milli¬ 
meter,  the  result  of  the  order  of  magnitude  of  the  currently 
assumed  values.  The  choice  of  0.1  mm  which  we  have  made  for 
the  £x  grid  would  hardly  permit  us  to  decrease  notably  the 
thickness  of  the  reaction  zone,  but  a  more  compact  grid  would 
assure  this  possibility. 

5.3.  CONCLUSION 

The  experimental  study  of  the  generation  of  detonation 
In  a  solid  granular  explosive  has  permitted  us  to  derive  a  few 
simplifying  assumptions  for  expressing  this  problem  in  the  form 
of  equations.  The  mathematical  model  prepared  in  this  way 
furnishes  results  which  are  in  rather  go^d  arreement  with  the 
experiments.  Nevertheless,  numerous  experiments  —  more 
accurate  than  those  described  here  --  and  more  complex  calcu¬ 
lations  will  be  necessary  to  illuminate  fully  this  delicate 
problem,  of  the  hydrodynamics  and  chemical  kinetics  of  explosive 
substances . 
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Diagrams  1 — 11  are  Interrelated  and  the  keys  are 
numbered  consecutively  as  follows: 

Key:  (l)  -  Generation  of  Detonation  by  Impact;  (2)  -  Profile 
of  the  pressure  signal  In  the  explosive  as  a  function  of  time; 
(3)  -  theoretical  results;  (A)  -  copper  projectile;  (5)  -  chem¬ 
ical  kinetics;  (6)  -  position  of  the  Impact  surface  at  the 
Initial  moment;  (7)  -  progress  diagram  of  the  reactive  shock 
front;  (6)  -  law  of  reaction  Inspired  by  holebumlng;  (9)  - 
free  surface;  f 1C)  -  solid  explosive;  (11)  -  law  of  reaction 
Inspired  by  grain-burning;  (12)  -  detonation  proiicto;  (13)  - 
reactive  shock  front;  (1A)  -  end  of  reaction;  (13)  -  retrograde 
detonation,  or  "retonation";  (16)  -  regardless  of  the  value 
of  p;  (17)  -  variation  of  the  reaction  rate;  (18)  -  Identical 
to  that  mentioned  on  diagram  7,  reference  i2];  (19)  -  effect 
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of  the  intensity  of  the  pressure  signal;  (20)  point  where  the 
velocity  discontinuity  is  situated;  (21)  -  initial  point  of 
the  establishment  of  the  steady  state  of  detonation;  (22)  - 
study  of  the  steady  state  of  detonation;  (23)  -  Von  Neumann's 
Feak  (theoretical);  (24)  -  mixture;  (25)  -  £as;  (26)  -  machine 
calculation;  (27)  -  theoretical;  (26)  -  dynamic  adiabatic  of 
the  solid  explosive  without  chemical  reaction  (machine  calcu¬ 
lation);  (29)  -  isentropic  originating  from  H  of  the  solid 
explosive  without  chemical  reaction  (machine  calculation). 
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